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ABSTRACT 


An experimental study of laminar flow heat transfer of an ethylene glycol/water 
mixture in an electrically heated horizontal tube using wire mesh (HEATEX) and 
twisted tape inserts was investigated. Twelve thermocouples, inserted in the tube wall 
at four longitudinal locations, enabled a mean inside experimental heat-transfer 
coefficient to be accurately measured. A constant wall heat flux boundary condition 
was placed on the wall by wrapping six 200 W flexible heater tapes tightly around the 
tube. The ethylene glycol/water mixture provided a coolant Reynolds number 
between 200-5000 and a Prandt! number between 30-140. Two smooth inside 
diameters and a roped tube profile were tested with and without the inserts. Heat- 
transfer correlations for tubes without inserts were ieatieae and compared with 
theory for both thermally and hydrodynamically developing flow. Correlations were 
also developed for the two types of inserts. Nusselt numbers for fully developed flow 
were found to be a function of Reynolds and Prandtl numbers for the wire mesh 
insert and a function of tape twist ratio, Reynolds and Prandtl numbers for the 
twisted tape insert. Heat transfer enhancements of — 7 for the wire mesh insert 
and over 4 for the twisted tape insert at high Reynolds numbers were obtained over 
the empty tube. By using these correlations in conjunction with earlier obtained 
refrigerant condensation data (using the same tubes, inserts, and coolant), more 


accurate values of the outside condensation heat-transfer coefficients were obtained. 
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NOMENCLATURE 


area, m? 

specific heat, J/kg K 

diameter, m 

Graetz number 

heat transfer coefficient, W/m*K. 
pitch for 180° rotation of tape 
current, amps 

thermal conductivity, W/m K 
length, m 

log mean temperature difference 
mass flow rate, kg/s 

Nusselt number 

perimeter, m or Power, W 
Prandtl number 

heat transfer rate, W 

heat flux, W/m? 

radius, m 

electrical resistance, ohms 


wall resistance, K/W 


XH 


Re Reynolds number 


Re, Reynolds number for swirl flow 
i Temperature, K 
t fin thickness, m 
U overall heat transfer coefficient, W/m? K 
u velocity component in x direction, m/s 
V velocity component in y direction, m/s 
V velocity, m/s or Voltage, volts 
x dimensionless axial coordinate, (x/D)/(Re Pr) 
y ° twist ratio, H/D 
: Greek Symbols 
5 tape thickness, m or boundary layer thickness 
E heat transfer enhancement ratio 
mn dynamic viscosity, kg/m s 
y kinematic viscosity, m?/s 
0 density, kg/m? 
Subscripts 
fe coolant or cross section 
D tube diameter | 
fd fully developed 
ee fin 
h hydrodynamically 


X1V 


sat 


t/c 


heated length 
inside or inlet 
laminar 

mean 

outside or outlet 
surface 
saturation 
thermally 
thermocouple 
local 


wall 





I. INTRODUCTION 


A. BACKGROUND 

With the increasing complexity of shipboard weapons and combat systems used 
aboard United States Naval vessels, the need for more energy efficient, lightweight, 
and high capacity refrigeration, air conditioning, and chilled water systems has grown 
tremendously. One way to meet this need is to redesign the refrigeration condenser. 
In redesigning these condensers, changes in condenser tube material can be 
combined with various ay transfer enhancement methods to achieve a lighter 
weight, more efficient mendenser tube. 

Emphasis on refrigerant condensation research reached a much higher level of 
importance due to the pressing environmental issues of ozone depletion and global 
warming. In September 1987, an international conference was held in Montreal, 
Canada, in which 24 nations representing the United Nation’s Environment Program 
| (UNEP), signed the Montreal Protocol on Substances that ee the Ozone Layer. 
This agreement called for a near term freeze on the production and consumption of 
chlorofluorocarbons (CFCs), the major contributor to ozone hi 
(Ref. 1]. In August 1988, the Environmental Protection Agency (EPA) 
adopted the regulations of the Montreal Protocol under the Clean Air Act [Ref. 2]. 


These regulations required a freeze in both production and consumption of CFC-11, 


12, 113, 114 and 115 at 1986 levels. The U.S. Navy currently uses CFC-12 and CFC- 
114 in its shipboard air conditioning and refrigeration systems. In June 1990 at a 
progress meeting held in London, England, UNEP delegates agreed to even larger 
reductions, calling for a 50% reduction by 1995, 85% reduction by 1997, and for 
CFCs to be totally phased out by the year 2000 [Ref. 3]. In the spring of 1992, 
President Bush moved up the CFC phase out deadline to 1995. 

Because of the Navy’s commitnient to the CFC phase out policy and the 
increasing need for lighter weight, more efficient refrigeration systems, refrigerant 


condensation research has surfaced to the forefront of Naval Research. 


B. REFRIGERANT CONDENSATION RESEARCH AT NPS 

In order to examine heat transfer characteristics of various condenser tube 
designs and alternative refrigerants, a condenser/evaporator test platform was 
constructed at the Naval Postgraduate School. Construction was begun by Zebrowski 
[Ref. 4] and completed by Mabrey [Ref. 5] in 1988. The test platform, shown in 
Figure 1.1, consisted of an evaporator and a condenser section with an associated 
ethylene glycol/water mixture coolant system. 

The evaporator section was designed so that in addition to supplying the 
refrigerant vapor to the condeneen it could also be utilized to conduct bundle boiling 
experiments. The condenser section consisted mainly of four instrumented condenser 
tubes and a series of auxiliary condenser coils used to control pressure in the 


apparatus. Details of design, construction and operation of the condenser/evaporator 


test platform are provided by Zebrowski [Ref. 4], Mabrey [Ref. 5] and Mazzone 
[Ref. 6]. 

Unlike a typical vapor compression refrigeration system used aboard ships, the 
test platform at NPS lacked a compressor between the evaporator and the condenser. 
This compressor is normally required to increase the pressure (and hence 
temperature) of the refrigerant vapor sufficiently enough so that it can give up its 
latent heat to the ambient sea water, the primary coolant. Because there was no 
compressor in the NPS test platform, the condenser must operate at the saturation 
conditions of the evaporator. Therefore, a much colder coolant, a refrigerated 
ethylene glycol/water mixture, was needed to condense the vapor, leading to laminar 
flow conditions within the condenser tubes. This in turn leads to poor heat transfer 
due to a large inside thermal resistance. : 

The primary objective of the condensation experiments of Mabrey (Ref. 5] and 
Mazzone [Ref. 6] was to accurately determine outside heat transfer coefficients. 
Briefly, this was accomplished by first determining the heat transfer rate to the 


* coolant by: 
g=mc, AT Cuil} 


where m is the coolant mass flow rate and AT the coolant temperature rise. 
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Figure 1.1 NPS Evaporator/Condenser Refrigeration Apparatus. 


The overall heat transfer coefficient, U,, could then be determined by: 


7 (1.2) 
A, LMTD 


where A, is the outside surface area of the tube and LMID is the log mean 


temperature difference, defined as: 


AT 


LMiD = 
een | (1.3) 
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The overall heat transfer coefficient can also be thought of in terms of a sum of 


thermal resistances. The overall thermal resistance is given by: 








ia ow (1.4) 
Ge Ay eA, 


The overall thermal resistance can be measured experimentally; the wall resistance, 
R,, is known, and the inside resistance can be determined from a chosen correlation. 


Therefore, the outside heat transfer coefficient, h,, can be calculated by: 


(1.5) 














There is one major drawback to this method. As mentioned earlier, when 
laminar flow heat transfer conditions exist within the condenser tube, it leads to poor 
heat transfer performance and therefore represents the dominant thermal resistance 


within the condenser. This inside. resistance then becomes the controlling factor in 


the determination of the outside heat transfer coefficient. Indeed, if the inside 
resistance is a large proportion of the overall resistance, then the accuracy of the 
calculated outside heat transfer coefficient becomes uncertain. Therefore, the 
accuracy of the outside heat transfer coefficient relies heavily on an accurate 
correlation to determine h;. 

An alternative is to reduce the inside controlling resistance by augmenting the 
inside heat-transfer coefficient. Augmentation techniques for lowering the inside 
thermal resistance primarily focus on repeatedly bringing fresh fluid from the coolant 
bulk stream to the heated inside surface. This results in breaking up the thermal 
boundary layer and thereby reducing the thermal resistance. The technique utilized 
in the NPS condensation experiments to break up this thermal boundary layer made 
use of two types of inserts, a twisted tape insert and a wire mesh (HEATEX) insert. 
These inserts and their use will be discussed in more detail in later chapters. 

The results of condensation experiments conducted by Mazzone [Ref. 6] show 
signs of inaccuracies related to the choice and accuracy of the inside heat transfer 
coefficient correlation. For example, in attempting to determine the outside heat 
transfer coefficient for a smooth tube with no insert, Mazzone found that the 
correlation used to determine the inside resistance yielded values that were greater 
than the actual measured overall resistance! Obviously, no sensible values of h, could 
be determined. Similar problems existed (but to a lesser extent, due to the greatly 
reduced inside thermal resistance) when using various types of inserts. Mazzone 


concluded that "... laminar flow correlations for twisted tape, HEATEX, and no 


insert conditions must be developed which have better accuracy than those used ... 
allowing for more accurate determination of outside heat transfer coefficients." 
Therefore, the determination of accurate inside heat transfer coefficients is 


paramount in the final outcome of the refrigeration condensation project. 


C. OBJECTIVES 
The main focus of this thesis was therefore to accurately determine inside heat 
transfer coefficients in support of refrigerant condensation experiments being 
conducted at NPS. 
The specific objectives of this study were to: 
1. Design and build a test apparatus to accurately determine inside heat transfer 
coefficients for the same tubes and inserts used in the condensation 
experiments. 


2. Develop correlations for inside heat transfer coefficients under laminar flow 
conditions both with and without inserts. 


3. Using the correlations mentioned above, reprocess existing condensation data 
to determine more accurate values of the outside heat transfer coefficients. 


II. LITERATURE SURVEY 


A. INTRODUCTION 
Laminar fully developed and developing forced convection flows have been 
analyzed in great detail for various boundary conditions. The study of heat transfer 
in laminar flow through a closed duct was first made by Graetz (Ref. 7] in 1883 and 
later by Nusselt [Ref. 8] in 1910. An excellent and thorough study of the theoretical 
development (with both analytical and experimental solutions) for various laminar 
forced convection internal flow boundary conditions is given by Shah and London 
(Ref. 9]. oo ae 
7 This survey will be concerned with laminar flow heat transfer in fully developed 
and developing flows subjected to both a constant wall temperature and constant 
heat flux boundary condition. It will also summarize some of the augmentation 
techniques used to enhance tubeside laminar flow heat transfer and discuss the 
details involved when using twisted tape and wire mesh inserts to enhance inside heat 


transfer. 


B. LAMINAR FLOW IN CIRCULAR DUCTS 
1. Boundary Layers and Entrance Lengths 


a. Hydrodynamic Boundary Layer 

For laminar flow in a circular tube, fluid enters the tube with a 
uniform velocity as shown in Figure 2.1; but when the fluid particles make contact 
with the tube wall, viscous effects cause them to assume a zero velocity. These 
viscous effects, associated with shear stresses acting in planes parallel to the tube 
wall, tend to retard the motion of particles in adjoining fluid layers, causing a 
momentum boundary layer to develop. The momentum boundary layer thickness, 6, 
increases until its growth is stopped by symmetry at the tube centerline. The viscous 
effects then extend over the entire cross section of the tube. The flow is then said to 
be hydrodynamically fully developed. At this point, the fluid velocity distribution at 
a given cross section is independent of the axial distance, x (i.e., du/dx=0) (Refs. 9, 
10, 11, 12}. | 

The hydrodynamic entry length for laminar flow (Re, s 2300) is 
defined as that region where the hydrodynamic boundary layer is developing and may 


be obtained from the approximation [Ref. 10]: 


= 0.05Re, Cc! 








This expression shows that for Rep=2300, the entrance length is 115 times the tube 
diameter, a considerable length for most experimental facilities. 


-Inviscid flow region 
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Figure 2.1 Laminar Hydrodynamic Boundary Layer Development in a Circular Tube 
(Courtesy [Ref. 10]) 


b. Thermal Boundary Layer 
When a fluid with a uniform temperature distribution enters a circular 
duct with a wall temperature that is greater than the uniform fluid temperature, 
temperature gradients in the fluid develop and convective heat transfer occurs, as shown 
in Figure 2.2. The region of fluid where these temperature gradients exist is called the 
thermal boundary layer. The thermal boundary layer thickness, 6, grows until a 
dimensionless fluid temperature distribution is independent of z The thermal boundary 


layer is then considered fully developed. This point is usually expressed in 


10 


Surface condition 





T(r, 0) Lp Tir,0) = T(r) 
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Figure 2.2 Thermal > en Layer Development in a Heated Circular Tube 
(Courtesy of [Ref. 10]) 


terms of a non-dimensional temperature profile given by [Ref. 11]: 


8 | 7)-TA)| _ (2.2) 
ex | T(x)-T_(x) 


where T, is the tube surface temperature, T is the local fluid temperature, and T, 
is the mean temperature of the fluid over the cross section of the tube. The exact 
shape of the fully developed temperature profile is dependent upon the boundary 
conditions imposed, either a uniform surface temperature or constant heat flux. 
For laminar flow, the thermal entry length may be expressed as the 


following [Ref. 10]: 


x | 
He = 0.05Re,Pr aa 


It should be noted that the rate of development of the velocity and 


temperature profiles in the entry region depend upon the fluid Prandtl number. For 


Pr = 1, both velocity and temperature profiles develop at the same rate. However, 
for Pr > 1, the hydrodynamic boundary layer develops more rapidly than the thermal 
boundary layer (X,3, < X,4,). For very large Prandtl] number (Pr > 100), X;,4, is so 
much smaller than X,,, that it is reasonable to assume that a fully developed velocity 


profile exists throughout the thermal entry region. [Refs. 9, 10, 11, 12] 


2. Fully Developed Flow 

Fully developed laminar flow of an incompressible fluid through a circular 
duct exists when both the thermal and velocity profiles are fully developed. Two 
boundary conditions must now be addressed, that of a uniform wall temperature and 
Rien wall heat flux. 

For these conditions, the general form of the energy equation in cylindrical 
coordinates for an incompressible fluid having constant properties (with no heat 
generation and neglecting body forces, pressure gradients and viscous dissipation), 
is given by (Ref. 11]: 


iid on eas (r] (2.4) 
ax or ror or}. 


It can be shown that for these two cases, the energy equation reduces to 


the following differential equations, 


(2.5) 





Constant Wall Heat Flux: i oO P oT _ il 
r or a 
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Uniform Wall Temperature: | 9 | oT ) malt, 
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For the constant wall heat flux boundary condition, Kays [Ref. 11] solves 
the differential equation (equation 2.5) for a fully developed velocity and 


temperature profile and expresses the Nusselt number as: 


oe be = 4.364 (2.7) 


Note that for this specified case, the heat transfer coefficient, h, depends 

only on the thermal conductivity, k, and the tube diameter, D. 
| A solution to the uniform wall temperature case is more complex, and 
involves evaluating the infinite series solution to the differential equation (equation 
2.6) for the hydrodynamically developed, thermally developing flow situation. In 


evaluating the series solution, the resulting asymptotic Nusselt number is given by 
lexef. 11]: 
Nu = 3.657 (2.3) 


This result is 16% less than the solution for the constant wall heat flux boundary 
condition. The asymptotic values in equations 2.7 and 2.8 are shown in Figure 2.3. 
Note that the inverse of the Graetz number is equivalent to the dimensionless axial 


distance x. 
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3. Entry Length Problem 

The solution to the energy equation for the entry region is far more 
difficult to obtain because velocity and temperature vary both axially and radially. 
There are two entry length problems to be considered, a combined entry length, 
where both temperature and velocity profiles develop simultaneously, and the 
thermal entry length problem, where the thermal conditions develop in the presence 
of a fully developed velocity profile. This latter case would exist if the place at which 
heat transfer begins is preceded by an unheated starting length, a situation that exists 
in most experimental facilities. It should be noted that for Pr << 1 (i.e., liquid 
metals), the thermal boundary layer would develop faster than the hydrodynamic 
boundary layer. However, this case has little practical application and is therefore 
rarely considered in the literature. 

Solutions have been obtained for both entry length conditions mentioned 
above [Refs. 9, 10, 11, 12] and these results are also presented in Figure 2.3. The 
Nusselt number is plotted against a dimensionless axial coordinate x (the inverse of 


the Graetz number) given by: 


x/D 


HF (2.9) 


Fully developed conditions are reached at x = 0.05 where the solutions for the 
thermal entry length converge with those for the combined entry length. For x" > 


0.05, the two types of entry length asymptotically reach the fully developed solutions 
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of 4.36 and 3.66 for the constant wall heat flux and uniform wall temperature 


boundary conditions respectively. 
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Figure 2.3 Local Nusselt Number Obtained from Entry Length Solutions for 
Laminar Flow in a Circular Tube (Courtesy [Ref. 10}) 


a. Thermal Entry Region - Uniform Wall Temperature 
The thermal entry length problem, first studied by Graetz in 1883 and 
later by Nusselt in 1910, is more commonly known as the Graetz Problem. They 
considered an incompressible fluid with constant properties flowing through a 
circular tube having a fully developed laminar velocity profile and a developing 
laminar temperature profile subjected to a uniform wall temperature. Shah and 
London [Ref. 9] present a numerically derived infinite series solution to the Graetiz 


Problem. 
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A more useful correlation presented by Incropera and Dewitt [Ref. 
10] for the uniform wall temperature thermal entry length problem was developed 
by Hausen and is of the form: 


0.0668(D/L)Re,,Pr 


Nu = 3.66 + ——______—_ 
1+0.04[(D/L)Re Phin 


(2.10) 


This correlation was used by Mazzone [Ref. 6] in his refrigerant condensation 
spaihaeds to predict the average inside heat transfer coefficient for a smooth tube 
without any type of insert. 

Shah and London [Ref. 9] present a more complete set of correlations 
derived from experimental data for heat exchanger tubes in laminar flow for local 
(Nu,) and mean (Nu,,) Nusselt numbers based on various x ranges. The correlations 


for local Nusselt number are given by: 


Nu, = 1.077(x«*)"“7-0.7 for x* = 0.01 (2.11) 
Nu_ = 3.657 +6.874(10%x°) Oe 7 for x* > 0.01 
and for mean Nusselt numbers by: 
Nu, = 1.615(x*)""3-0.7 for x* < 0.005 
Nu =1.615(x%*)"-0.2 for 0.00S<x* < 0.03 (2.12) 
Nive eoes memes for x* = 0.03 
= 


The reader is reminded that both the Hausen correlation and the 
Shah and London correlations presented above are for a thermally developing flow 


with a fully developed velocity profile subjected to a uniform wall temperature. 
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b. Thermal Entry Region - Constant Wall Heat Flux 
For the thermal entry length problem under a constant wall heat flux 
boundary condition, Shah and London [Ref. 9] present several numerical solutions 


but recommend the following approximate equations for the local Nusselt number: 


Nu, = 1.302(x*)"-1 fomeon sox 1057 
Nu, = 1.302(x*)"-05 for 5x10% <sx° < 15x10 — (2.13) 
Nu_ = 4.364+8.68(10%x °)%e 4 (ome os 105 


Shah and London also recommend the following approximate equations for the mean 


Nusselt number: 


Nu, = 1.953") for x* = 0.03 
2.14 
Nu, = 4.364+20722 for x* > 0.03 iets 





The above two sets of equations are for a thermally developing, 


hydrodynamically developed flow subjected to a constant wall heat flux. 


c. Combined Entry Length 
| As discussed earlier, the rate at which the temperature and velocity 
profiles develop is dependent upon the Prandtl iia At medium to high neu 
numbers (Pr > 100), the velocity profile develops much more rapidly than the 
temperature profile and a fully developed velocity profile assumption is valid; for 
such conditions the equations given above should be used. 
For 1 < Pr < 100, however, a fully developed velocity profile is no 


longer a valid assumption and a combined entry length solution must be developed. 
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Again, Shah and London [Ref. 9] present a comprehensive study of various solutions 
to the combined entry length problem, including two solutions for local Nusselt 
numbers developed by Churchill and Ozoe [Ref 13], the first, for a uniform wall 


temperature: 


Nu, +1.7 
5.357{1 +[(388/17)x°]PP* 
fk (2:15) 
43 
ar/(284x°) 
{1 +(Pr/0.0468)*7}"7{1 +[(388/2r)x *} PP" 


Bytes 





and the second, for a constant wall heat flux: 


Nu_+1 
5.364{1 +[(220/ar}x°] “PR 
; 3710 (2.16) 


$3 
q/(115.2x°) 
{1 +(Pr/0.0207)*9}'7{1 +[(220/ar)x *}-1P 


ie: 





From the preceding paragraphs, it can be seen that the subject of 
laminar flow forced convection in circular ducts has been studied in great detail 
yielding both analytical and experimental results. The reader is referred to Shah and 
London [Ref. 9] and Kays and Crawford [Ref. 11] for a more complete and thorough 


review of this subject. 
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C. AUGMENTATION TECHNIQUES TO ENHANCE TUBESIDE HEAT 
TRANSFER 


1. General Introduction 

As discussed earlier, when laminar flow heat transfer occurs within tubes, 
it usually represents the dominant thermal resistance in a heat exchanger. These 
large resistances result in low inside heat transfer coefficients. In recent years, the 
requirement for more efficient heat exchangers has stimulated considerable interest 
in techniques to augment or enhance tubeside heat transfer. 

The various methods to enhance tubeside heat transfer can be classified 
into two categories, passive and active techniques. Passive techniques involve no 
external energy (except increased pump work) to produce the necessary 
enhancement. Passive techniques are divided into the following methods: 

1. Internally finned tubes (longitudinal or spiral fins) designed to increase the 
inside heat transfer surface area. 


2. Surface roughness techniques aimed at agitating flow rather than increasing 
heat transfer surface area. 


3. Swirl flow techniques which involve the use of inlet vortex generators, 
periodically spread propellers or twisted tape inserts designed to cause 


enhancement by creating a rotating and/or secondary flow. 


4. Displaced promoters which alter the fluid flow near the surface by the use of 
inserts such as wire mesh, static mixer elements, rings or disks. 


5. Coolant additives designed to increase the conductivity of the coolant. 


6. Compound techniques that involve more than one of the above mentioned 
methods. 
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Active techniques utilize an external energy source to promote the 
required augmentation. These techniques include: (1) mechanical aids, (2) heated 
surface vibration, (3) fluid pulsation, (4) electrostatic fields, and (5) suction or 
injection. 

All of the above techniques are aimed at creating turbulence in the 
relatively slow moving or stagnant boundary layer fluid by displacing this fluid from 
the heated surface and muang it with the bulk fluid. exert in all of these 
techniques is a balance between increased heat transfer rate and increased pumping 
power requirements/costs. Bergles and Joshi [Ref. 14] present an excellent review 
of augmentation techniques and compare heat transfer data taken at low Reynolds 
number flow subjected to both uniform watt temperature and constant wall heat flux. 

Of particular interest to this thesis are two passive techniques, the swirl 
flow technique using twisted tape inserts and the displaced promoter technique using 
wire mesh inserts. These two methods of enhancement of tubeside heat transfer are 


discussed in more detail in the following sections. 


2. Heat Transfer Enhancement with Twisted Tape Elements 
Among the most common swirl flow augmentation techniques is the use 
of twisted tape inserts. The tapes usually consist of a thin stainless steel, brass or 
copper strip which is the width of the inside diameter of the tube. The tapes are 
twisted by clamping one end to an overhead, attaching weights to the other end and 


then twisting to the desired specifications [Ref. 15]. Enhancement occurs not only 
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due to the increased path length of flow, but also due to secondary flow effects and 
fin effects caused by tape contact with the inside tube wall. 

Analytical studies of twisted tape inserts by Date and Singham (Ref. 16] 
suggest that swirl flows improve heat transfer by as much as a factor of 70; however, 
no experimental data have confimned their unrealistically large predictions. The first 
experimental work of swirl flow using twisted tapes was conducted by Hong and 
Bergles (Ref. 15]. Their primary objective was to develop experimentally based 
correlations for predicting heat transfer coefficients for laminar flow of water and 
ethylene glycol in an electrically heated tube (heated length of 1.22 m, inside 
diameter of 10.2 mm) with two sizes of twisted tape inserts. They showed that the 
Nusselt number was independent of axial location, suggesting that fully developed 


conditions existed. Therefore, the mean Nusselt number should take the form: 
Nu = f(Re,,Pr,y). (Zale) 


Re, is the Reynolds number for the swirl flow given by: 
Renee ac aa (2.18) 
7 (D,-45) 


and 6 is the tape thickness. The twist ratio, y, is defined as: 


H | 
= Si De 19 
y D (Zao) 


as 


where H is the pitch for a 180° twist of the tape. Hong and Bergles then combined 
Re, and y as a single parameter and came up with the following correlation for 


values of Re,/y > 10: 


Naz = US83 Pr (heyy) aie (2.20) 


For values of Re,/y < 10, the following correlation was obtained: 


Nu, = 5.172{1 +5.484 x 10° Pr°"(Re /y)!9)5 (2.21) 


It should be noted that all fluid properties in equations 2.20 and 2.21 
should be evaluated at the bulk fluid temperature. One additional important 
observation was made by Hong and Bergles regarding tube wall temperature. They 
found a significant variation of temperature at one axial location with two different 
tape orientations at that position of 90° and 120°, suggesting that the temperature 
profile was related to the tape orientation. Hong and Bergles concluded that for a 
fully developed flow, inside heat transfer coefficients depend on twist ratio and could 
be improved by a factor of two or three over an empty tube. These conclusions were 
confirmed by later studies of Saha et. al. [Ref. 17] and Marner and Bergles (Ref. 
18]. 

3. Heat Transfer Enhancement with Wire Mesh (HEATEX) Elements 

The use of wire mesh inserts to promote inside heat transfer falls into the 

displaced promoter technique mentioned earlier. One such wire mesh insert was 


commercially developed by CAL-GAVIN Limited (UK) and called "HEATEX." 
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This HEATEX insert consists of a central wire core on which a series of 
wire loops or petals are attached. The petals are attached to the core at an inclined 
angle that face the oncoming fluid flow. Ideally, each petal makes contact with the 
tube wall in the form of an arc, thereby maximizing disturbance of fluid along the 
tube wall. In addition to the boundary layer disturbance, radial mbang occurs caused 
by repeated redistribution and ming of fluid through the mesh of loops near the 
tube axis [Ref. 19]. Various wire loop densities and nominal diameters are 
available, depending on the application. 

Initial studies by Gough et. al. (Ref. 19] with HEATEX inserts suggest 
enhancements of a factor greater than ten times that of a similar tube with no insert. 
Additional experiments were conducted by Oliver and Aldington [Ref. 20] in which 


a correlation was determined: 


Nu = 0.232Re 4 Pr °-* Ces), 


It must be remembered, however, that this correlation is for a specific loop density 
and nominal diameter HEATEX element. 

The two types of HEATEX mixing elements supplied by CAL-GAVIN 
(nominal diameters of 13.26 and 10.16 mm) did not match the loop densities of 
previous studies; therefore new correlations needed to be developed. Mazzone [Ref. 
6] developed correlations for these inserts from limited data supplied by CAL- 
GAVIN. For the 13.26 mm and 10.16 mm elements, these correlations were 


respectively: 


Nu = 0.226Re°SPr° (13.26mm) (2.23) 


and 


Nu = 0.063Re°*Pr* (10.16mm) (2.24) 


Since these correlations were developed from a very limited amount of 
data taken under conditions that did not exactly match that of Mazzone’s 
condensation experiments, they were always subject to imaccuracy. Therefore, the 
need for more accurate correlations under the exact conditions used in the 


condensation experiments are required. 
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Il. EXPERIMENTAL APPARATUS 


A. SYSTEM OVERVIEW 

The test apparatus utilized is shown in Figures 3.1 and 3.2. The apparatus is a 
closed loop, low pressure system utilizing an ethylene glycol/water mixture as the 
coolant. The coolant is circulated by a centrifugal pump from a supply tank, through 
a reducing manifold and flowmeter into the test section. Coolant is then returned to 
the supply tank via a return header. 

The test section consisted of an inlet and outlet thermocouple well/mvang 
chamber to measure inlet and outlet coolant temperatures and the test tube itself. 
The tube was instrumented with twelve thermocouples (embedded in the tube wall) 
to determine inside wall temperatures. A constant heat flux boundary condition was 
achieved using a series of flexible electrical heating elements wrapped around the 
tube. Monitoring the system as well as data acquisition and data processing were | 


provided by a Hewlett-Packard microcomputer and Data Acquisition System. 
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Figure 3.1 Primary and Secondary Coolant Systems. 
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Figure 3.2 Data Acquisition and Control Systems. 
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B. TUBES AND INSERTS TESTED 


1. Tubes Tested 

Three types of tubes were utilized in the expenments, a copper smooth 
tube, a copper/nickel externally finned tube, and a copper/nickel corrugated or roped 
tube. These tubes correspond to the tubes utilized in the refrigerant condensation 
experiments of Mabrey [Ref. 5] and Mazzone [Ref. 6]. 

The smooth copper tube provided a baseline set of experimental data 
which were then compared with previous experimental and theoretical work for the 
inside heat transfer coefficient conducted by Shah and London [Ref. 9] for no insert, 
Oliver and Aldington (Ref. 20] for the HEATEX insert, and Hong and Bergles [Ref. 
15] for the twisted tape insert. 

The Cu/Ni finned tube provided a smaller inside diameter to be used in 
these comparisons. Note that the fins play little part in the experiments. The wall 
thermocouples were buried in the tube wall and were used to measure the inside wall 
temperature, which had a smooth bore. Consequently, for the purposes of these 
experiments, the tube behaved like a smooth copper tube with a different internal 
diameter. 

The final tube tested was a Cu/Ni corrugated (or roped) ane. 
commercially referred to as a "Korodense” tube. The tube was manufactured and 
provided by Wolverine Tube Company of Decatur, Al. All tube dimensions and 


specifications were taken from the Wolverine Technical Bulletin No. 4020 [Ref. 
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21]. Figure 3.3 is a schematic of the Korodense tube. Tube specifications for all 


three tubes are listed in Table E.1 of Appendix E. 
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Figure 3.3 Wolverine Korodense Tube (Courtesy [Ref. 21]) 


Zz. Insert Elements Tested 

Two types of insert mvang elements were utilized in this investigation, a 
twisted tape insert and a wire mesh insert commercially referred to as a "HEATEX" 
tube insert. In both cases, two sizes were required due to the two different tube 
inside diameters used. 

The twisted tape inserts were wail of a strip of brass 0.559 mm thick and 
either 10.16 mm or 13.26 mm wide. The tapes were manufactured in the NPS 
Mechanical Parncenae machine shop by clamping them to the overhead, attaching 
a weight to the opposite end, and twisting them to the desired twist ratio. A twist 
ratio of 3 and 4 corresponding to the two respective tube inside diameters were 
manufactured. Once installed in the tube, the edges of the tape should ideally have 


been in contact with the tube wall. 
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The second type of insert used was a commercially available product 
referred to as HEATEX radial mixing elements. These elements were supplied by 
Cal Gavin of Birmingham, U.K. They are manufactured from stainless steel and 
consist of a central wire core onto which a series of wire loops or petals are attached. 
Each petal was inclined at an angle facing the oncoming flow. The outside diameter 
of the elements were slightly oversized for the tube, ensuring that close contact with 
the tube wall in the form of an arc, rather than a point, was achieved. As with the 
twisted tape elements, different diameter HEATEX elements were needed, 13.26 
mm used in the smooth and Korodense tubes and the other 10.15 mm for the Cu/Ni 


finned tube. Figure 3.4 shows both the twisted tape and HEATEX inserts. 
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Figure 3.4 Twisted Tape and HEATEX Inserts 
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C. COOLANT SYSTEM 


1. Primary Coolant System 

An ethylene glycol/water mixture (54/46% by weight) was stored in a 1.5 
m?’ insulated supply tank. The mixture was circulated through the primary system by 
a 1% hp, constant speed, centrifugal pump. The coolant was discharged from the 
pump via a 76 mm diameter PVC pipe to a manifold where it was reduced to a 15.9 
mm diameter Tygon flexible tube. A rotameter type flowmeter with a ball valve at 
its entrance was used to control the flow rate through the system. Flowmeter 
' calibration results are discussed in Appendix A. 

Flow entered the test section through a 90° bend which housed the 
thermocouple well and exited through a mixing chamber/thermocouple well 
combination (see Figure 3.1). Prior to entering the exit mixing chamber, the coolant 
passed through a 10 cm length of flexible tubing containing a high density wire mesh. 
This wire mesh, in conjunction with the mixing chamber, served to minimize radial 
temperature Carations and provide a uniform outlet coolant temperature. The flow 
then exited the test section through more flexible tubing to a manifold and return 


PVC header, suspended above the apparatus, and back to the coolant supply tank. 


2. Secondary Coolant System 
In order to achieve a wide range of Reynolds number, coolant inlet 
temperatures were varied from -10°C to +20°C. This was accomplished by a 


secondary coolant system comprised of a 28 kW (8 ton) external refrigeration system. 


ep 


A 30 gpm, .75 hp centrifugal pump continuously erculaeed coolant between the 
supply tank and a chiller barrel where it was cooled to the desired temperature. 
Temperature was maintained by a thermostatically controlled solenoid valve at the 
discharge side of the chiller barrel. This system allowed for a range of coolant 
temperatures between -20°C and ambient to be selected while maintaining the 


coolant inlet temperature approximately constant (+2°C) during each run. 


D. INSTRUMENTATION AND TEST SECTION 


1. Instrumentation 

All temperature measurements were made using copper/constantan type-T 
teflon coated thermocouples and were previously calibrated by Mabrey (Ref. 5]. The 
inlet coolant temperature was measured by a single thermocouple while the outlet 
coolant temperature was measured by two thermocouples located immediately after 
the mvang chamber. | 

Each test tube was fitted with twelve wall thermocouples to determine the 
average inside wall temperature. Thermocouples were placed at four different 
longitudinal positions along the axis of the tube. Three thermocouples at each 
location were spread evenly around the tube at 120°C intervals. Figure 3.5 shows a 
typical tube with the exact longitudinal and nodal thermocouple locations. Actual 
longitudinal locations varied slightly for each tube and are given in Table E.1 as x/D, 
where x is the length from the start of the heated length to the thermocouple 


location and D, in the inside diameter of the tube. The two extreme longitudinal 
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locations of the thermocouples corresponded to the beginning and end of the 
"effective" condensing length of the tubes being used in the condensation experiments 


of Mabrey [Ref. 5] and Mazzone (Ref. 6]. 


INSTRUMENTED TUBE 
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Figure 3.5 Instrumented Tube 
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2. Thermocouple Attachment Method 

In order to attach the thermocouples to the wall of the tube, short 
longitudinal grooves were cut at each location. These longitudinal grooves were cut 
to a depth of % the tube wall thickness. It was decided that a thermocouple 
embedded within the tube wall vice being attached to the surface would lead to a 
more accurate tube wall temperature and minimize the effect of direct contact 
heating from the heating elements. Thermocouples were then placed in tnese grooves 
and either soldered in place or tack welded and covered with a liquid metal epoxy. 
Different methods for attaching and covering the thermocouples were used for each 
tube based on the tube material and configuration titended Below 

In the case of the smooth copper tube, the thermocouples were-silver 
soldered in place and the remainder of the groove then covered with silver solder. 
Since the tubes were copper, this was relatively easy to accomplish. Care was taken 
to ensure that the thermocouple remained in Pp with the tube wall and that no 
air pockets were created (which would cause an additional large resistance to heat 
flow). Also, great care was taken to prevent any melting of the teflon insulation past 
the thermocouple joint. This condition could have caused an electrical contact 
between the two wires of the thermocouple prior to the tube wall, leading to an 
unrealistically high observed wall temperature. 

A different method was required to attach the thermocouples to the 
copper/nickel finned tube due not only to the external fins but also the 10% nickel 


content. Initially, a small diameter hole was machined into the tube to minimize fin 
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damage (rather than cutting a longitudinal groove). Since the thermocouple was to 
be embedded half way into the tube wall, this hole proved to be relatively deep. 
Also, because of the nickel content of these tubes, the preheat temperature required 
to ensure that the silver solder would adhere to the tube wall was much higher. 
These higher required temperatures and the hole depth caused a significant amount 
of the thermocouple teflon insulation to be melted, allowing for potential contact 
between the thermocouple wires at some point other than the tube wall. This 
technique was therefore abandoned in favor of the longitudinal groove method. 
However, thermocouple insulation was again being destroyed due to the higher 
required preheat temperatures. An alternative solder, soft solder, with lower preheat 
and melting temperatures, was therefore used. The final combination of longitudinal 
grooves aaa soft solder was successful. 

For the copper/nickel Korodense tube, the groove and soft solder method 
was again tried. Initial tests in the apparatus showed wall temperatures at one 
longitudinal position varying by as much as 15°C. On closer examination of the 
thermocouple connections (using a magnifying glass), it was revealed that insulation 
damage was evident on over half of the thermocouples. Several more attempts at the 
groove and soft solder method were tried, all leading to the same results. 

A third method nae a tack welding machine (acquired and adapted 
for our use) was tried. Thermocouples were tack welded into the bottom of the 
groove resulting in no damage to the insulation and a liquid steel epoxy was used to 


fill the remainder of the groove, providing additional strength and protection to the 
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connection. This final method proved to be the most successful and is recommended 


for use in any follow-up experiments. 


3. Test Section 

In order to achieve a constant heat flux, each tube was tightly wrapped 
with six silicone rubber flexible electrical heating elements rated at 240 V, 200 W 
each. Each strip was 2.5 cm wide and 50.6 cm long. The heating elements were 
connected in parallel to a 240 Vac power supply via a plug panel and power 
switchboard. Power to the six heaters was controlled by a STACO, 240 Vac, 23.5 kVa 
rheostat (see Figure 3.2). 

To minimize heat loss from the heater to ambient, the tube was wrapped 
in a 20 mm thick neoprene type fnsutWOH: This type of insulation proved to be 
extremely effective giving a negligible heat loss. It also provided insulation against 
heat influx from ambient at the colder coolant inlet temperatures. Actual heat loss 
and heat influx experiments and their results are discussed in more detail in 
Appendices B and C. All other tubing and PVC piping was wrapped in similar type 
insulation. | 

The whole test section was then placed in a restraining stand which 
minimized outside vibration effects and ensured that each tube remained horizontal 
and in the same relative position (to the remainder of the apparatus) for each 


experiment. 
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E. DATA ACQUISITION AND CONTROL SYSTEM 

A Hewlett-Packard 9125 microcomputer (MC) was used to control a Hewlett- 
Packard 3497A Automatic Data Acquisition Unit (DAU) (see Figure 3.2). The DAU 
read the output from the twelve tube wall thermocouples and three coolant 
thermocouples (one inlet and two outlet). Additional channels were assigned to sense 
heater input voltage and current at the plug panel. The voltage was converted to 
heater power using the heater resistance and the calibration discussed in Appendix 
D. 

Thermocouple readings were made in microvolts which were converted to 
temperatures (along with the heat flux calculations) using the data reduction 
program. The data reduction program, DRPSING, will be discussed in more detail 
in a later chapter. All temperatures, power, and heat flux measurements were 
monitored using the MC and DAU. DAU channel assignments are listed in Table 


Bel. 
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TABLE 3.1 DATA ACQUISITION UNIT CHANNEL ASSIGNMENTS 


_ CHANNELNO. =| ME NO. "MEASURED QUANTITY 


Thermocouple Pos 1 
Thermocouple Pos 2 
Thermocouple Pos 3 
Thermocouple Pos 4 
Thermocouple Pos 5 
Thermocouple Pos 6 
Thermocouple Pos 7 
Thermocouple Pos 8 
Thermocouple Pos 9 
Thermocouple Pos 10 
Thermocouple Pos 11 
Thermocouple Pos 12 
Tube Inlet 
Tube Outlet 
Input Current 
Input Voltage 
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IV. EXPERIMENTAL PROCEDURES 


A. CALIBRATION EXPERIMENTS 

Several calibration experiments were required prior to performing the main 
experiments. Figure 4.1 depicts a typical tube being tested. The primary heat source 
is from the heating elements (Gunn). However, to complete the heat balance on the 
system, several other heat sources and heat sinks needed to be accounted for. The 
first experiment conducted was to determine the amount of heat being lost from the 
heating elements to the atmosphere vice going to the coolant. Heat loss experiments 
are discussed in Appendix B. 

The second experiment was to account for two additional potential heat 
sources: heat leakage into the coolant from ambient (4q,,,,) and frictional temperature 
rise (q;,,.)- Approximately 8 cm at either end of the tube, along with the flexible 
tubing section prior to the thermocouple wells, were not covered by the heating 
elements. Since coolant inlet mmpertanres were primarily below ambient 
temperature, a potential for heat leakage to the coolant from ambient existed. Also, 
because of these cold inlet temperatures, frictional heating could exist, especially 
when using the various inserts. Heat leakage into the system and frictional effects 


were considered together as one heat influx. These experiments were conducted with 


the heater tapes turned off; any temperature rise was then monitored using the inlet 
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and outlet thermocouples. Appendix C discusses the experiments conducted to 
determine this part of the heat balance equation. As mentioned earlier, the test 


section and all other tubing were insulated to minumize these effects. 


HEAT BALANCE 





Figure 4.1 Heat Balance. 


B. -~ PROCEDURES FOR DETERMINATION OF INSIDE HEAT TRANSFER 
RATE 7 


1. General Procedures 

Each tube utilized in this study was tested with no insert and the 
appropriate twisted tape and HEATEX inserts corresponding to the tube inside 
diameter. For each insert condition, ‘aut nominal coolant inlet temperatures, -10°C, 
0°C, +10°C and +20°C were chosen to increase the range of experimental Reval 
and Prandit numbers that could be attained. To check that heat flux variations did 
not affect the results, data were taken at heat fluxes of 10 kW/m? and 15 kW/m? for 
each inlet temperature. For each insert, coolant inlet temperature and heat flux 
setting, coolant velocity was varied between 0.2 m/s and 1.4 m/s. Since the inside 


diameter of the Cu/Ni finned tube was smaller than the other two tubes, higher 
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coolant velocities of up to 2.4 m/s could be attained. A matrix of the individual data 
runs and data file names is shown in Table 4.1. 

With the tube installed and all calibration experiments completed, the 
experiments to determine the inside heat transfer coefficient were conducted. Since 
flowmeter settings for specific coolant velocities were based upon the inlet coolant 
temperature, the system was initially operated at a 100% flowmeter setting with no 
power to the heating elements for 15 minutes. This allowed the inlet thermocouple 
well to stabilize. Once the desired inlet coolant temperature was reached and 
maintained, power to the heating elements was set to the predetermined heat flux 
setting for that data run (10 or 15 kW/m’). The flowmeter was then adjusted to 
achieve the desired coolant velocity based upon-flowmeter calibration results 
(discussed in Appendix A). Flow was maintained for approximately 10 minutes to 


ensure steady state equilibmum conditions prior to taking data. 


TABLE 4.1 DATA FILE NAMES 


TEMP (C)/ 
cree ataet Tape Insert 
Tse fomoe| = [moe ese [ene ay [| re [om 
a a ad ee yt tl 
+15 +15 +15 
| surtor_| CFTTO2 |_cettog | crnxoe | xonzo2 | Ree | 
ae ne 
| oy10_|_swntoe | _surroz_| sunxo2 | crnros | cerros | cenxos | xonros | korres | KoHxos | 
| ois __| sunzoz |_surroz | swnxoz | centoa | cerros | cewxos | xowzos |_xorros_| xowxog_ 
| -sosro_| sunzor |_suttoa_| swaxos | cento2 | cerroa | crnxoa | xonto | xotroa | KoHxo4 
|_-10/18_{ swwtor | _swrroa | sunxoa | cenroz | certoa | crnxoa | xontoa | xortos | xonxoa | 
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Velocities were increased from 0.2 m/s to 1.4 m/s in steps of 0.2 m/s and 
then decreased from 1.3 m/s to 0.3 m/s in similar steps. This method of increasing 
on even values and decreasing on odd values was done not only to determine if any 
hysteresis effects existed (i.e., a difference in the increasing and decreasing data), but 
also to ultimately increase the number of data points acquired. Also, since data were 
immediately plotted (as Nu, versus Re), the method allowed for an easy 
determination as to whether a data point had reached a steady state condition. 

Inlet and outlet coolant temperatures, tube wall temperatures, heat flux 
and heating element power settings were al] continuously monitored using the 
computer program DRPSING in conjunction with the Data Acquisition System 
discussed earlier (DRPSING is listed in Appendix G). Individual data runs for each 


tube and insert condition were conducted in a similar manner. 
2. Specific Experiments 


a. Copper Smooth Tube 
The first tube tested was the copper smooth tube. This tube was 
chosen to provide a baseline set of data needed for validation of system performance 
and for comparison with data acquired using other tubes. It should also be noted that 
two complete independent sets of data covering all insert conditions and inlet 
temperatures were taken over a period of two weeks. Data from both sets were 
compared and found to be in very good agreement (+2%); this validated the 


repeatability of data obtained from the system under the same conditions. 
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The initial intention for all tubes was to start with a given insert at 
the highest coolant inlet temperature, take data at each of the heat flux settings and 
then proceed to the next coolant temperature. Once one insert was completed, the 
system could then be shut down and the insert changed. This procedure would then 
be repeated for all three insert conditions. However, the coolant supply tank and 
coolant system were used concurrently for other experiments which required the 
coolant temperature to be -20°C. Since it took approximately 36-48 hours for the 
coolant temperature to increase from -10°C to +20°C (the sump was not fitted with 
a heater), only the copper smooth tube was done in this way. The method for the 


other two tubes is discussed below. 


b. Cu/Ni Externally Finned Tube 

The second tube tested was the Cu/Ni finned tube. For a given inlet 
temperature and tube insert, the same procedure as with the copper smooth tube was 
used. However, rather than decreasing coolant temperature with the same insert, the 
coolant temperature was maintained and the insert changed until all three of the 
insert conditions were completed for a given inlet temperature. Once all the data for 
a given inlet temperature had been taken, the supply tank temperature was 
decreased to the next inlet temperature condition. It took only about one hour for. 
the secondary coolant system to decrease the temperature by 10°C. Since it was a 
relatively short task to change inserts, the 36 to 48 hour waiting period required to 


raise coolant temperatures back to +20°C between insert changes was eliminated. 
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It should be noted that for the Cu/Ni finned tube, the highest nominal 
inlet temperature attained was only +15°C vice +20°C. Colder than normal ambient 
temperatures existed preventing supply tank temperatures from reaching +20°C. 
Also, due to the smaller inside diameter, coolant velocities could be varied between 


0.3 m/s and 2.4 m/s. 


c. Cu/Ni Korodense Tube 

The same procedure used with the finned tube was also cede the 
Korodense tube. In addition, another test was conducted with this tube. It was noted — 
that larger than normal wall temperature vanations existed at one longitudinal 
location when using the two inserts. It was suspected that the point of contact 
between the insert and tube wall was directly below a thermocouple position causing 
a higher wall temperature at that point. To verify this, the system was shut down by 
securing the coolant pump and the heating elements so that the flowmeter and 
power settings remained exactly the same. The outlet end of the tube was opened 
and the insert moved downstream by approximately 2 cm. The system was then 
closed and the coolant flow and power to the heating elements restored. Wall 
temperature variations still existed but were significantly less. This procedure was 


repeated several times with similar results. 


V. DATA REDUCTION 


A. AVERAGE INSIDE HEAT TRANSFER COEFFICIENT 

The main objective of the data reduction scheme was to determine the average 
inside heat transfer coefficient (h,). This in tum was used to determine a mean 
Nusselt number (for the whole heated length of the tube) to be used as the basis for 
establishing the desired correlations. 

To accomplish this task, the Data Acquisition Unit discussed in Chapter III 
made use of the computer program DRPSING. A listing of the program is given in 
Appendix E. The program allowed the user te monitor information from the 
instrumentation attached to the apparatus, acquire and then process the data to 
determine the desired output. In doing this, the user would select the tube and insert 
condition as well as the nominal! coolant inlet temperature for that particular data 
run. The desired coolant velocity was then selected. Based on the actual measured 
coolant inlet temperature, coolant mixture properties and flowmeter calibration 
results (see Appendix A), the appropriate flowmeter setting (as a %) was calculated 
and set on the flowmeter. 

The program then monitored coolant inlet and outlet temperatures, the twelve 
wall thermocouple temperatures and the actual measured heat flux and power 


settings. Once all readings reached a steady state condition (approximately 10 min.), 
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the user could terminate the display session and acquire data for that coolant velocity 
setting. 

The DAU measures emf from the thermocouples and then converts these 
values to temperatures using the temperature calibration equation developed by 
Mabrey [Ref. 5]. Heater power was determined by measuring the voltage at the plug 
panel and calculating the power based on heater resistance and the power calibration 
contained in Appendix D. 

The twelve wall thermocouples (placed axially and radially) were used to 
minimize the uncertainty in determining the average inside wall temperature. These 
twelve wall temperatures were averaged and corrected (to account for the depth of 
burial using a simple conduction calculation) to determine the average inside wall 
temperature, T,,. 

The physical and thermodynamic properties of the ethylene glycol/water 
mixture (density, dynamic viscosity, thermal conductivity, heat capacity) were 
determined using equations and figures from Cragoe [Ref. 22] and Gallant (Ref. 
23] and correlated by Mabrey [Ref. 5]. From these quantities, the coolant mass 


flow rate was calculated by: 
m = AeA (5.1) 


where A, is the inside cross sectional area of the tube and V, is the set coolant 


velocity. The coolant Reynolds number, Re, was calculated by: 
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V. Dip ¢ 
H. 


Re (5.2) 


D 


where D,; is the inside diameter of the tube. The dimensionless distance (x) in the 


flow direction is calculated as: 


D. 
ie rae (5.3) 
Re, Pr 


where x is the axial distance from the beginning of the heated length. 

To determine the heat transfer rate to the coolant, several factors were 
considered. The basic equation to determine the heat transfer rate is given by: 

q = mc,AT (5.4) 
where AT is the measured coolant temperature rise along the tube. The temperature 
rise required is that for the actual heated length of the tube that is subjected to the 
constant heat flux conditron, AT,,,. To determine AT,,, from AT, corrections have to 
be made to account for heat loss, heat leakage and heat transferred to the unheated 
ends. Figure 5.1 depicts a typical tube and the correction terms considered in the 
heat balance. 

The overall AT is first corrected for heat loss to the ambient, (q,,,,) and heat 
leakage into the system, (4,,,,). The heat leakage term consisted of any temperature 
rise due to the influx of heat from ambient as well as that caused by frictional 
effects. Appendix D discusses the experiments conducted to determine both this heat 


loss and heat leakage correction. 
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TEMPERATURE RISE CORRECTIONS 


T 
t/c.cut 
@ 





Figure 5.1 Temperature Rise Corrections. 


To take account of the heat transferred to the coolant from the unheated ends, 
the heat transfer rate to the coolant over the entire length of the tube (heated and 
unheated) is first calculated. The average heat transfer coefficient can then be 
determined by: | 

h, = aes (55) 
where A, is the inside surface area based on the heated length (L,,) and LMTD is 
the log mean temperature difference defined as: 

At this point in the calculation, however, the heat transfer rate to the coolant 


from the heated length alone (q,,,) is not known. Since the heating elements did not 
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(9 a ala i 
T,,-T,, (5.6) 
In 
7 -T ug 


cover the entire length of the tube, account must be taken of heat added to the 
coolant from the unheated lengths at the inlet and outlet of the tube. This was 
accomplished by considering these ends as annular fins. An adiabatic fin tip and 
outer wall boundary condition was applied and the heat transfer from the fins, q, (as 


given by Incopera and Dewitt [Ref. 10]), determined from: 


q, = hi PRA, @,tanhmL, 


where L, = fin length 
P = fin perimeter = wDi 


A, = fin cross sectional area 


7 2 Zz 
Seale § Da D) 


9, = T,-T., (for inlet fin) 


= Jue ou. (fOr outlet fin) 


h 0.5 
m= iP 
kA 


¢ 
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The fin calculation requires a known average inside heat transfer coefficient (h,). 
Therefore, an initial estimate of h, is calculated by: 

hi. = ais (5.8) 
where q,,,. is known. The amount of heat that passes through the unheated ends is 
then calculated using the fin correction above. The heat transfer rate to the coolant 


from the heated lengtu (q,,,) is then calculated by: 


Gis = dea (5.9) 
A new average inside heat transfer coefficient (h,) can then be calculated using 
equation 5 5 and an iteration process carried out until convergence to a final average 
inside heat transfer coefficient is reached. 

An additional correction for the Cu/Ni finned tube is needed to account for the 
change in inside diameter of the smooth ends (needed to fit through the tube sheet 
in the condenser). The heat transfer coefficient based on the smaller inside diameter 
corresponding to the finned portion of the tube is what ultimately is desited As 
shown in Figure 5.2,-the heated length extended beyond the smaller diameter section. 
To correct for this, the annular fin correction discussed earlier was applied to 
' determine the temperatures at the inlet and outlet of the heated length. The coolant 
temperature rise along the heated length was assumed to be linear, and therefore the 
temperature rise along the smooth larger inside diameter lengths (L, and L,) could 


be determined by the simple ratio: 
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sa rt (5.10) 





eee and Al = 


HL HL 


where L,, is the heated length of the tube and ATj,q, the temperature nse along the 
heated length. These corrections were then added or subtracted from the corrected 


finned inlet and outlet temperatures. 


Cu/Ni EXTERNALLY FINNED TUBE 
TEMPERATURE RISE CORRECTIONS 


INLET oa ae lee OUTLET 


Figure 5.2 Cu/Ni Finned Tube Temperature Rise Corrections. 
With the average inside heat transfer coefficient known, the mean Nusselt 


number can be determined by: 
ee: (5.11) 


where k is the thermal conductivity of the ethylene glycol/water mixture. 
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The experimentally determined mean Nusselt number was then compared to 
various theoretical and experimental correlations, based on the insert condition. 
These correlations were discussed earlier but are briefly mentioned here. 

| For the no insert condition, the data were compared to two theoretical 


correlations, the Hausen correlation [Ref. 10] given by: 


NR Pate si tin (5.12) 
1+0.04[(D/L)Re, Pr] 


and the Shah and London correlation [Ref. 9] given by: 


Nu = 1.93(x*)4 (5.13) 
The reader is reminded that the Hausen correlation is for a uniform wall 
temperature boundary conditron and was used in Mazzone’s condensation 
experiments [Ref. 6]. The Shah and London correlation is for a constant wall heat 
flux condition. It should also be noted that equation 5.13 is used since x’ for the test 
section was always less than 0.03. 
For the twisted tape insert the experimental mean Nusselt number was 


compared to the correlation developed by Hong and Bergles (Ref. 15]: 


Nu, = 0.383Pr°>5(Re /y)° (3.14) 


Equation 5.14 is used because values for Re,/y were always greater than 10. 


Sy 


For the HEATEX insert, the data were compared to the two correlations 
developed by Mazzone [Ref. 6]. For the 13.26 mm inside diameter tube, the equation 


is given by: 
Numesn0226Ke0 Fr -— (5.16) 


and for the 10.16 mm inside diameter tube by: 


Nu = 0.063Re°”Pr° (5.16) 


B. ADDITIONAL TUBE LENGTHS 


In addition to the above calculations for the entire tube length, mean Nusselt 





numbers for three additional tube "lengths” could be determined. This was 
accomplished by sectioning the heated length of the tube into four cells as shown in 
Figure 5.3. Each cell contained one longitudinal thermocouple position. Three 
additional tube lengths were then determined by adding cell lengths together. The 
first length, L,, consisted of cell 1 alone, L, consisted of cells 1 and 2 and L, 
consisted of cells 1, 2 and 3. The fourth length was ee the entire heated length 
as previously discussed. 

The coolant temperature rise along the heated length of the tube was assumed 
to be linear. The correction to the wall heat flux to account for variation in electrical 
resistance with temperatures was assumed to be negligible. Therefore, the 


temperature rise along each cell could be determined by: 
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SECTIONED TUBE 


be 
T/C POS 1 T/C POS 2. T/C POS 3 T/C POS 4 


— — 


Figure 5.3 Additional Tube Lengths. 












(5.17) 





Average inside wall temperatures were based on the thermocouples within that tube 
"length." Therefore, the. temperature rise, heat flux and average inside wall 
temperatures were all known for each tube length. Applying the same calculation 
method as before, the mean Nusselt number for these three additional tube “lengths” 
could be calculated. 

This adaption to the program 2 not done until after the smooth tube 
experiments had been completed; however, this extra data were taken for the 


remaining two tubes. 
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C. LOCAL NUSSELT NUMBERS 

The program was also adapted to determine local inside heat transfer 
coefficients at the four axial thermocouple positions. Corrections discussed earlier 
to determine the temperatures at the inlet and outlet of the heated length were first 
accounted for. Exact thermocouple longitudinal positions were entered into the 
program as L/D ratios so that a local x could also be calculated. Again, the 
_ temperature rise along the heated length was assumed to be linear; therefore, the 
bulk coolant temperature at the particular thermocouple position could be 
determined by the simple ratio method discussed earlier (equation 5.10). 

The local heat transfer coefficient, h,, was then determined by: 


h q” 5.13 
gt (5.13) 





where ii. is the average inside wall temperature based on the three thermocouples 
at that position and T, is the bulk coolant temperature at that position. The local 


Nusselt number was then given by: 


Nu = eam (5.19) 





The local Nusselt number for all four thermocouple positions was then determined. 
For the smooth tube experiments, local Nusselt numbers were calculated only 
for position 1. The program was later adapted to include local calculations at 


positions 2, 3, and 4 for the Cu/Ni finned tube and Korodense tube. 


VI. RESULTS AND DISCUSSION 


A. GENERAL 

The experimental heat transfer results are initially presented in terms of the 
mean Nusselt and Reynolds numbers for all insert conditions. From these graphs, it 
can be seen how the Nusselt number is inlet temperature dependent, indicating the 
need for a Prandtl number (raised to an exponent) in the final correlation. However, 
although not specifically shown, the data shows no dependence on the heat flux 
setting. In addition, the Nu vs. Re graphs for the no insert condition show when a 
" transition to turbulent flow occurs. 

The heat transfer results are then compared to existing correlations based upon 
the insert condition. For the no insert condition, the Nusselt number is plotted 
against x (dimensionless x coordinate). Note that for a given tube, L/D is constant 
and x’ varies due to Re and Pr changes. The data are then compared to the Hausen 
(Ref. 10] correlation (equation 2.10) and the Shah and London [Ref. 9] correlation 
(equation 2.14). Recall that the Hausen correlation is for a uniform wall temperature 
boundary condition and was used by Mazzone [Ref. 6] in the refrigerant 
condensation experiments. Since the present data were taken under a constant heat 
flux boundary condition, comparison with the Shah and London correlation was also 


required. The twisted tape insert data are compared to the correlation developed by 
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Hong and Bergles [Ref. 15] (equation 2.20) which was also used by Mazzone. The 
HEATEX insert data are compared to the two correlations developed from the 
limited data received from CAL-GAVIN by Mazzone (equations 2.27 and 2.28). 
Correlations are developed for each tube and insert condition based upon a 
| least squares linear regression. An r-squared value, a measure of the accuracy of the 
least squares curve fitting to the data, is used to determine the uncertainty of the 
correlation. The correlations and the corresponding r-squared values are summarized 
in Table 6.1. Note that the 0.2 m/s and 0.3 m/s data were eliminated prior to 
developing all correlations due to the high degree of uncertainty in these data points. 


TABLE 6.1 SUMMARY OF CORRELATIONS 


~ ~~ SUMMARY OF CORRELATIONS 


NO INSERT TWISTED HEATEX 
| TAPE 


Nu = Nu = Nu = 
oo 4.1387Xstar*-0.234 .512Pr~.35(Res/y)*.588 | .184Pr*.46Re*.641 


r-sq = 71.9% r-sq = 99.5% r-sq = 98.9% 


Nu * Nu = Nu = 
1.84Xstar~-0.341 .148Pr*.35(Res/y)”.687 | .126Pr~.46Re~.657 


r-sq = 93.0% r-sq = 96.7% r-sq = 95.9 


Nu = Nu = Nu # 
2.047 Xstar°-0.328 | .589Pr°.35(Res/y)*.498 | .282Pr°.46Re~.499 


r-sq 7 89.2% r-sq 2 98.9% r-3q = 94.3% 
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B. SMOOTH TUBE 


1. Smooth Tube - No Insert 

The first tube tested was the copper smooth tube with no insert. Heat 
transfer data is presented in Figure 6.1 in terms of the Nusselt number as a function 
of Reynolds number for each nominal inlet temperature. The first thing to notice is 
the sudden increase in Nu at a Re number of about 2300. This is indicative of a 
change from laminar to transition flow where the heat transfer is greatly increased 
due to improved mixing. This contradicts Mazzone’s [Ref. 6] assumption of a 
transition to turbulence occurring at Re > 4000. The data in the transition region 
was primarily for the +20°C and +10°C nominal inlet temperatures. Figure 6.1 also 
shows distinct groupings of data at given inlet temperatures, indicating a Prandtl 
number effect. These groupings are shown more clearly in Figure 6.2 which represent 
data under laminar flow conditions only (Re s 2300) and Figure 6.3 for data in the 
transition region only (Re > 2300). | 

To develop a useful heat transfer correlation, the mean Nusselt number 


is plotted as a function of x in Figure 6.4. Recall that x" can be determined by: 


xXx° = L/D _ (2:2) 
Re, Pr 


For a given tube, L/D is a fixed quantity. Therefore, for a given Prandtl number 


(essentially fixed inlet temperature), x is dependent only upon the Reynolds number. 
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Figure 6.4 also shows a comparison of the experimental data with the 
Hausen correlation (equation 2.10) for a uniform wall temperature (used by 
Mazzone) and the Shah and London correlation (equation 2.14) for a constant heat 
flux. The Hausen correlation appears to significantly underestimate the heat transfer 
coefficient found in the present experimental data, i.e., the present experimental data 
indicate a lower inside thermal resistance. This is consistent with Mazzone’s [Ref. 6] 
results in that the Hausen correlation was predicting larger inside resistances than 
the overall measured resistance. However, it must always be remembered that the 
Hausen correlation is for a uniform wall temperature boundary condition and this 
Manidlindeed lead to lower heat-transfer coefficients. 

The data also fall slightly above (but parallel to) that predicted by the 
Shah and London correlation. Due to the constant heat flux boundary condition used 
by Shah and London, one may expect better agreement (as seen). However, the fact 
that the data fall slightly above can be attributed to secondary flow induced by the 
taratus inlet conditions (90° bend due to the thermocouple well). Also, as before, 
a departure from laminar flow can be Pegi: Figure 6.4. Since the transition data 
forms two distinct lines (+20°C and +10°C), it indicates that x is not a good 
correlating: parameter for transitional flow. 

After removal of the transitional flow data, a correlation was developed 


based upon the form of the Shah and London correlation. The resulting present 
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correlation for the mean Nusselt number for laminar flow with no insert is given by: 


Nae = 4137 aye (6.1) 


Equation 6.1 and the Hausen correlation, equation 2.10, are both shown 
in Figure 6.5 with the present laminar flow data. The general form of Figure 6.5 is 
explained by the fact that at low x (high Re) the heat transfer is improved due to 
slightly better mixing. As x increases (coolant velocity decreases), the heat transfer 
decreases. 

Figure 6.5 also shows a high degree of scatter in the data. Initially, it was 
thought that this scatter was due to the flow not being hydrodynamically fully 
developed. However, it is seen to be more pronounced at large x. Recall that for a 
given nominal inlet temperature, large x represents a low Reynolds number. For a 
low Reynolds number, the hydrodynamic entry length (given by X/D = 0.05Re) is 
the shortest. Also, since we are dealing with relatively large Prandtl numbers (30- 
140), the assumption that the velocity profile is fully developed throughout the 
thermal entry region should be fairly valid. Therefore, it is felt that this scatter is not 
due to the flow being hydrodynamically wndeveloped! 

This scatter is believed to be due to the effects of the different inlet 
temperatures, indicating that x’ is not an ideal correlating parameter. That is to say, 
the Prandtl! number exponent needs to be adjusted to account for the different inlet 


temperatures. It should also be noted that due to the high degree of detail in Figure 


60 


6.5, the scatter appears more pronounced. However, in Figure 6.4 with larger scale 


increments, the scatter seems minimal. 


2. Smooth Tube - Twisted Tape 
Results for the smooth tube with a twisted tape insert are shown in 
Figures 6.6 and 6.7. A definite Pratidt] number dependence can be seen in Figure 6.6 
where the mean Nusselt number is presented as a function of Reynolds number. 
Following the guidance of Hong and Bergles (Ref. 15] the dimensionless parameters 
Re, and y were combined as one parameter Re,/y and the mean Nusselt number and 


0.35 


Prandtl number were combined as Nu/Pr~>. The final correlation shown in Figure 


6.7 is: 


Nu = 0.512Pr°5 (Re jy)“ (6.2) 


This correlation is in close agreement with the Hong and Bergles correlation 
(equation 2.20). The slight differences can be attributed to differences in the 
experimental apparatus (including the insert itself) and the possible existence of 
secondary flows. Figure 6.7 shows little inlet temperature dependence, indicating that 
the Prandtl number exponent of 0.35 is probably fairly accurate. Furthermore, there 
is little experimental data scatter, giving increased confidence in this correlation. 
3. Smooth Tube - HEATEX Insert 

HEATEX insert results for a smooth tube are shown in Figures 6.8 and 

6.9. Again, on a Nu vs. Re basis, it can be seen from Figure 6.8 that there is a 


definite dependence on the inlet temperature. Using the same form as the 
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correlation of Oliver and Aldington [Ref. 20], the resulting correlation for the mean 


Nusselt number is: 


Nu = 0.164Pr°“Reo™! (6.3) 


This is compared in Figure 6.9 to the correlation developed and used by Mazzone 
(equation 2.23). It is clearly evident that Mazzone’s correlation underpredicts the 
inside thermal resistance and therefore gives a larger inside heat transfer coefficient. 
Again, Figure 6.9 shows little dependence on inlet enierature indicating that the 
Prandtl number exponent is fairly accurate. Comparisons of the enhancementin heat 
transfer gained by the use of HEATEX and twisted tape insert, over the no insert 


condition are given section F of this chapter. 
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Figure 6.1 Smooth tube - No insert (Nu, vs. Re) results. 
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Figure 6.2 Smooth tube - No insert (Nu,, vs. Re) laminar flow results. 
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Figure 6.3 Smooth tube - No insert (Nu,, vs. Re) transitional flow results. 
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Figure 6.4 Smooth tube - No insert (Nu,, vs. x) results. 


66 


SMOOTH TUBE - NO INSERT 


LAMINAR FLOW 


Mean Nusselt 


Nu = 4.14Xstar*-0.234 





0.0005 0.001 0.0015 0.002 0.0025 0.003 0.0035 
Xstar 


, t20C ,+10C , OC . -10C 


Nu vs Xstar 
SMNIXNULAM 


Figure 6.5 Smooth tube - No insert laminar flow correlation results. 
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SMOOTH TUBE - TWISTED TAPE 


Nu vs Re 
os 
a2 
e r 
o a 
o 
“ a 
Ps a3 
A 
a ad .3 
Ne Pd “ 
a® o 
a as 
a. ? 
& 
ve ¢ ’ 
. a 
p w 
. — ee 
n % 
4s 3 ~ 
8 
4 8 


500 1000 1500 2000 2500 3000 3500 
Reynolds 
@ a ] 0C & OC A -|] OC 


Figure 6.6 Smooth tube - Twisted tape (Nu,, vs. Re) results. 


68 
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Figure 6.7 Smooth tube - Twisted tape correlation results. 
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SMOOTH TUBE - HEATEX 
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Figure 6.8 Smooth tube - HEATEX (Nuzg, vs. Re) results. 
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Figure 6.9 Smooth tube - HEATEX correlation results. 
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C. CU/NI EXTERNALLY FINNED TUBE 


1. Cu/Ni Finned Tube - No Insert 

The second tube tested was the Cu/Ni externally finned tube. Recall that 
this tube has a smaller inside smooth diameter than the smooth tube (10.15 mm vs. 
13.26 mm). The no insert results for the mean Nusselt number as a function of 
Reynolds number are shown in Figure 6.10 and in more detail in Figures 6.11 and 
6.12 for laminar and transitional flows, respectively. It can be seen that a slightly 
higher Reynolds number (Re = 2700) for the departure from laminar flow exists for 
this tube. Again, the dependence on inlet temperature is clearly evident. 

The data are compared to the Hausen correlation (equation 2.10) and the 
Shah and London correlation (equation 2.14) in Figure 6.13. As expected, the 
Hausen correlation again underpredicts the inside heat transfer coefficient. However, 
the laminar flow data is in closer agreement with the Shah and London correlation 
than they were for the smooth tube. This can be attributed to the different diameter 
affecting the thermal boundary layer development. The resulting correlation for the 


laminar flow mean Nusselt number is shown in Figure 6.14 and is given by: 


Nu, = 1.84(x*)-0 (6.4) 


The same conclusions about data scatter can be drawn as for the smooth tube 


results, although in general the scatter seems to be less. 


wz 


2. Cu/Ni Finned Tube - Twisted Tape Insert 
Figures 6.15 and 6.16 show the results for the Cu/Ni finned tube using a 
twisted tape to enhance tubeside heat transfer. The inlet temperature dependence 
in Figure 6.15 seems to be reduced, but still evident. The form of the correlation was 
again that of Hong and Bergles’ correlation. It was developed in the same manner 


as for the smooth tube and is given by: 


Nu, = 0.148Pr°"* (Re jy)?’ (6.5) 


-The data and resulting correlation are compared to the Hong and Bergles correlation 
(equation 2.20) in Figure 6.16 where the latter is seen to overpredict the mean 
Nusselt number. This difference can be attributed to several factors. A small 
diameter tube with a larger diameter smooth entrance length could cause severe 
secondary flows. Also, because of the smaller diameter tube a different size twisted 
tape (twist ratio, y=3) was used, which had a slightly looser fit than the tape used 
_ with the smooth tube. Additionally, Prandtl number effects need to be investigated 


more closely, as Figure 6.16 shows a definite trend with inlet temperature. 


3. Cu/Ni Finned Tube - HEATEX Insert 
Data for the Cu/Ni finned tube using the HEATEX insert are presented 
in Figures 6.17 and 6.18. Figure 6.17 (Nu vs. Re) again shows a dependence of the 
data on the coolant inlet temperature. Since the HEATEX element used for this 
tube was a smaller diameter (10.16 mm) with a different loop density than that used 


for the smooth tube, a new correlation needed to be developed. Again, using the 
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form of the correlation proposed by Oliver and Aldington, the final proposed correlation 


is given by: 
Nu, = 0.126 Pr° Re (6.6) 


The data and resulting correlation are compared to the correlation developed 
by Mazzone (equation 2.24) for this insert in Figure 6.18. It can be seen that Mazzone’s 
correlation is adequate for predicting the inside heat transfer coefficients at low Reynolds 
number (Re < 1400), where both correlations merge. However, at higher Reynolds 
numbers, Mazzone’s correlation overpredicts the present experimental data. The 
differences in the two correlations are primarily due to the limited amount of data 
available to Mazzone when developing his correlation. Furthermore, the daia he used 
were taken on a different experimental facility (by the manufacturers). As are before, 
the data shown in Figure 6.18 still indicate a dependence on inlet temperature and hence 
the value of 0.46 for the Prandtl number exponent needs to be reinvestigated and 


optimized. 
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Figure 6.10 Cu/Ni externally finned tube - No insert (Nu,, vs. Re) results. 
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Figure 6.11 Cu/Ni externally finned tube - No insert (Nu,, vs. Re) laminar flow 
results. 
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Figure 6.12 Cu/Ni externally finned tube - No insert (Nu,, vs. Re) transitional flow 
results. 
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Figure 6.13 Cu/Ni externally finned tube - No insert (Nu,, vs. x) results. 
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Figure 6.14 Cu/Ni externally finned tube - No insert laminar flow correlation results. 
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Cu/Ni FINNED TUBE - TWISTED TAPE 
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Figure 6.15 Cu/Ni externally finned tube - Twisted tape (Nu,, vs.- Re) results. 
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Figure 6.16 Cu/Ni externally finned tube - Twisted tape correlation results. 
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Cu/Ni FINNED TUBE - HEATEX 
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Figure 6.17 Cu/Ni externally finned tube - HEATEX (Nu, vs. Re) results. 
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Figure 6.18 Cu/Ni externally finned tube - HEATEX correlation results. 
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D. KORODENSE TUBE 


1. Korodense Tube - No Insert 

The final tube tested was the Korodense tube. This tube was first tested 
with no insert. The results for the mean Nusselt number as a function of Reynolds 
number are shown in Figure 6.19 and in more detail for laminar and transitional flow 
in Figures 6.20 and 6.21, respectively. As was seen in the smooth tube case, a 
transition to turbulent flow occurred at a Reynolds number of approximately 2300. 

Figure 6.22 compares the data to Hausen’s correlation (equation 2.10) and 
the Shah and London correlation (equation 2.14). As before, the transitional flow 
data are seen as two distinct lines for data taken at +20°C and +10°C nominal inlet 
temperatures. Under laminar flow conditions, the data are seen to fall very closely 
to Shah and Londons’ correlation, but significantly above Hausen’s correlation. This 
was also found for the other two tubes. This result was expected since a mildly 
corrugated wall is not good in enhancing laminar flow and the tube behaves 
essentially like a smooth tube. Significant enhancement for Korodense tubes is not 
expected until Rep > 10° (Ref. 21]. 

The resulting correlation for laminar flow data is shown on Figure 6.23 


and given by: 


Nae 205 (6.7) 


The Hausen correlation (equation 2.10) used by Mazzone again under predicts the 


present measured inside heat transfer coefficient. It was expected that the mean 
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Nusselt number would be the same (or maybe slightly higher) as for the smooth tube 
case for the reasons already mentioned. However, the Korodense data fall slightly 
below the smooth tube data and closer to Shah and Londons’ smooth tube 
correlation (equation 2.14). No explanation can be given for this result. Data scatter 


explanations given for the other tubes also apply for the Korodense tube. 


2. Korodense Tube - Twisted Tape Insert 

The results for the Korodense tube with a twisted tape insert are shown 
in Figures 6.24 and 6.25. The same twisted tape that was used for the smooth tube 
was used for the Korodense tube. Note that the Korodense tube nominal inside 
diameter (13.35 mm) is slightly larger than the smooth tube inside diameter (13.26 
mm). This small difference in inside diameters (together with corrugation in the 
Korodense tube) led to a looser fit for the twisted tape insert. 

The mean Nusselt number as a function of Reynolds number is shown in 
Figure 6.24 which again clearly indicates the inlet temperature dependence. The data 
were again correlated using the form of Hong and Bergles’ correlation. The present 
data were compared to their correlation (equation 2.20) in Figure 6.25. The 


correlation developed for the present data is: 


Nu, = 0.569Pr°* (Re Jy)°*8 (6:8) 


This correlations falls below the Hong and Bergles correlation (equation 6.20) as well 
as for the smooth tube. This result is thought to be primarily due to the looser fit of 


the twisted tape inside the Korodense tube. The corrugations also leave areas of no 
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contact between the tape and the inside wall of the tube. Saenant flow or even air 
bubbles may form in these gaps, especially at low Reynolds numbers, causing an 
interruption of the swirled flow. The thermal boundary layer, therefore, is not as 
efficiently mixed as in the smooth tube. For these reasons, lower inside heat transfer 


coefficients than those assumed by Mazzone should be expected. 


3. Korodense Tube - HEATEX Insert 
The last tube-insert condition vested was the Korodense tube with the 
HEATEX insert. The same HEATEX insert used for the smooth tube was used 
here. Figure 6.26 shows the mean Nusselt number as a function of Reynolds number 
and again shows an inlet temperature dependence. The resulting correlation for the 


present data ts given by: 


Nu, = 0.282Pr°“Re%™ (6.9) 


The data and this correlation are compared with the correlation Mazzone used 
(equation 2.23) in Figure 6.27. It is evident that Mazzone’s chosen correlation 
overpredicts the inside heat transfer coefficient Erion This again is thought 
to be primarily due to the incomplete mixing of the thermal boundary layer caused 


by non-contact regions existing between the tube wall and HEATEX petals. 
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Figure 6.19 Korodense tube - No insert (Nu,, vs. Re) results. 
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Figure 6.20 Korodense tube - No insert (Nu,, vs. Re) laminar flow results. 
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Figure 6.21 Korodense tube - No insert (Nu,, vs. Re) transitional flow results. 
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Figure 6.22 Korodense tube - No insert (Nu, vs. x ) results. 
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Figure 6.23 Korodense tube - No insert laminar flow correlation results. 
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KORODENSE TUBE - TWISTED TAPE 
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Figure 6.24 Korodense tube - Twisted tape (Nu, vs. Re) results. 
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Figure 6.25 Korodense tube - Twisted tape correlation results. 
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Nu = 0.569Pr*.35(Res/y)*0.498 
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Figure 6.26 Korodense tube - HEATEX (Nu, vs. Re) results. 
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KORODENSE TUBE - HEATEX 
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Figure 6.27 Korodense tube - HEATEX correlation results. 
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E. INSERT COMPARISON BETWEEN TUBES 

Figures 6.28, 6.29 and 6.30 each compare the three correlations for each tube 
for no insert, twisted tape, and HEATEX inserts, respectively. It can be seen that 
with no insert (Figure 6.28), the Korodense correlation is about 10% below the 
smooth tube data. Although one may expect no increase in the heat transfer with a 
roped tube under laminar flow conditions, one certainly would not expect to find 
such a decrease with a similar diameter tube. This difference, therefore, can only be 
attributed to experimental uncertainty, probably in the wall temperature 
measurement. The finned tube correlation is lower than both the Korodense (=5%) 
and smooth tube (= 15%) correlations; this tube had a smaller inside diameter and 
it may be that different induced flow patterns are set up enue a reduction in heat 
transfer. Experimental uncertainty must also not be ruled out as to the cause of the 
difference. 

Figure 6.29 shows the same three tube correlations, but for the twisted tape 
insert. The smooth tube again has the highest performance. The finned tube is =20% 
below at a high Re number and nearly 40% lower at a low Re number. This could 
be due to a different twist ratio used for the finned tube due to the smaller inside 
diameter. The Korodense tube is similar to the finned tube, giving slightly lower heat 
transfer at high Re number and slightly higher heat transfer at low Re. The reason 
for this tube being significantly lower than the smooth tube is though to be due to 


the tape not contacting the inner surface of the tube as well as for the smooth tube, 
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thereby leaving regions of stagnant flow. These regions would tend to inhibit mixing 
and lead to a lower heat transfer performance. 

For the HEATEX comparison (Figure 6.30), similar conclusions can be made, 
although the finned tube seems to have better relative performance indicating that 
the small diameter HEATEX insert used in the small diameter tube is more effective 
than the twisted tape. The difference between the smooth and Korodense tubes are 


of similar magnitude to that shown in Figure 6.29 (about 40% at high Re). 
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Figure 6.28 No insert comparison under laminar flow for all three tubes. 
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Figure 6.29 Twisted tape insert comparison under laminar flow for all three tubes. 
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Figure 6.30 HEATEX insert comparisons under laminar flow for all three tubes. 
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F. LAMINAR FLOW ENHANCEMENT RATIO 

The comparisons made in Figures 6.28 through 6.30 are made more 
complicated by the fact that the uncertainty associated with the wall temperature 
measurement is unknown and most likely different for each tube. For this reason, it 
is not known how much of the divergence seen in these figures is due to differences 
in thermocouple attachment (giving uncertainty to the wall temperature 
measurement) and how much is an actual flow induced phenomenon. A comparison 
of the relative effects of the three inserts for a given tube is more accurate, however, 
due to the fact that any uncertainty for a given tube wall thermocouple attachment 
combination would be the same. 

To determine the degree of tubeside heat transfer enhancement achieved by 
the use of the twisted tape and HEATEX inserts under laminar flow, an 
enhancement ratio of the mean Nusselt number for the insert condition (Nu,,,.,,) to 
that for the no insert condition (Nu,, ..,.,) Was calculated. The enhancement ratio is 


therefore defined as: 


Figures 6.31 through 6.33 present mean Nusselt number as a function of Reynolds 
number under laminar flow conditions for each tube under the three insert 
conditions. For the purposes of this study, enhancement ratios were calculated at the 


two extreme Reynolds numbers of 200 and 2300 with a constant Prandtl number of 
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70 using the correlations developed in the previous sections. Table 6.2 summarizes the 
enhancement ratios for each tube and insert condition. 

It is evident from these figures that the greatest enhancement is achieved using the 
HEATEX insert. An enhancement factor of over seven was achieved for the smooth tube 
using the HEATEX insert at Re = 2300. This reduced to about 2.6 at low Re. 
Significantly smaller enhancements were achieved for the twisted tape inset for the 
smooth tube, indicating that for laminar flow, the stripping and mixing effect of the 
HEATEX element is considerably more effective than the simple swirl flow effect 
provided by the twisted tape. It should be mentioned here that HEATEX was developed 
specifically with laminar flow enhancement in mind. In turbulent flow (Re > 10°), the 
enhancements obtained are much smaller (around 1.5 to 2) and both twisted tape and 


HEATEX elements give similar enhancement values. 


TABLE 6.22 SUMMARY OF LAMINAR FLOW ENHANCEMENT RATIOS 
(defined by equation 6.10) 











SMOOTH TUBE FINNED TUBE KORODENSE TUBE 
TWISTED TWISTED TWISTED 
TAPE HEATEX TAPE HEATEX TAPE HEATEX 
Low Re 
(Re = 200) 
High Re 
(Re = 2300) 7.19 6.71 2.93 4.25 
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For the finned and Korodense tubes, the HEATEX element also provides the 
highest enhancements. At high Re, the enhancements provided by the twisted tape 
and HEATEX inserts are lower than those obtained with the smooth tube for the 
reasons mentioned previously. However, at low Re, the enhancements for the three 
tubes are about the same for a given insert, indicating the "stagnant" regions 


mentioned earlier have less of an effect on the total heat transfer. 
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Figure 6.31 Insert condition comparison for laminar flow in a smooth tube. 
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Figure 6.32 Insert condition comparison for laminar flow in a Cu/Ni externally 
finned tube. 
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Figure 6.33 Insert condition comparison for laminar flow in a Korodense tube. 
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G. REPROCESSED REFRIGERANT CONDENSATION DATA 

Figure 6.34 shows the condensation data of Mazzone [Ref. 6] with the outside heat 
transfer coefficient recalculated using the present correlations. The data presented are for 
the top copper smooth tube. Similar recalculations can also be carried out for the other 
two tubes (Cu/Ni finned and Korodense tubes) as well as for condensation bundle data. 
The data of Mazzone are reprocessed using the Modified Wilson Plot Technique (see 
Mazzone [Ref. 6] for details). This technique essentially lets the leading coefficient of 
the correlation ’float’ to try to smooth out any uncertainty in the correlations being used. 
As previously stated, the chosen correlation for the no insert condition (even when 
reprocessing using the Wilson Plot Technique) gave values of inside thermal resistances 
greater than the measured overall thermal resistance, thereby yielding negative values of 
the outside heat transfer coefficient, h°. This explains why Figure 6.34 has no data for 
Mazzone’s no insert condition. 

It can be seen that by using the present developed correlations for the no insert 
condition, the calculated values of h, are significantly improved (i.e. positive values). 
In addition, for the other two inserts, the new correlations increase h, by about 10% for 
the twisted tape insert and nearly 20% for the HEATEX insert. The newly calculated 
values of h, when using these inserts are now very close and just above the Nusselt 
prediction as may be expected for low vapor velocity film condensation. Mazzone 


commented that his data were below Nusselt by about 10% 
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and thought that this was a consequence of the length of his condenser tubes violating the 
uniform wall temperature assumption of the Nusselt theory. The fact that these data are 
so close gives confidence in the new correlations, as the insert itself should have no 
influence on the magnitude of h,. By significantly reducing the inside resistance, the 
insert simply increase the accuracy of h,. Of course, the no insert condition should also 
give similar values of h,. Figure 6.34 shows that they are too high, indicative of the fact 
that the correlation for the no insert condition, although much improved, is still not ideal. 
If one recalls the scatter obtained in the present data for the no insert case, together with 
the fact that the inside resistance is still dominating, then the results shown in Figure 
6.34 are not too surprising and still very encouraging. It is recommended, therefore, that 
all future refrigerant condensation experiments carried out on the bundle facility should 
be conducted with the inserts using the newly formulated correlations. 

One final note is the fact that the present correlations were formulated 
using a constant wall heat flux, whereas for the condensation experiments, one 
should ideally use correlations developed under uniform wall temperature conditions. 
One may expect that the inside heat transfer coefficient would increase slightly under 
constant wall heat flux conditions, thereby decreasing the outside heat transfer 
coefficient. However, from Figure 6.34 it is clear that the use of the constant wall 
heat flux correlations seems to increase h,. Thus it appears that with the inserts, the 


inside resistance is so reduced that the boundary condition imposed has little 
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influence. Future investigators should look into this boundary value problem in 


more detail. 
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Figure 6.34 Reprocessed refrigerant condensation data using new correlations. 
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VII. CONCLUSIONS 


The following conclusions are based upon results obtained in this study. 
A. GENERAL 


i. A test apparatus was successfully designed and built to accurately determine 
inside heat transfer coefficients for the same tubes and insert conditions used 
in the condensation experiments. 


2. Correlations to accurately determine inside heat transfer coefficients under 
laminar flow conditions both with and without inserts were successfully 
developed. 


3. Reprocessing of Mazzone’s [Ref. 6] refrigerant condensation data using these 
new correlations has yielded more accurate values of the outside heat transfer 


coefficient. (Note: the name of the program used for the reprocessing is 
DRPCON?7.) 


B. SPECIFIC 


1. A departure from laminar flow into a transitional flow exists at Reynolds 
numbers of approximately 2300 for the smooth and Korodense tubes and 2700 
for the Cu/Ni finned tube. This contradicts Mazzone’s [Ref. 6] assumption of 
laminar flow up to a Reynolds number of approximately 4000. 


2. Laminar flow results for the no insert condition were in close agreement with 
the Shah and London [Ref. 9] correlation for a constant wall heat flux 
boundary condition for all three tubes tested. 


3. The Hausen [Ref. 10] correlation used by Mazzone [Ref. 6] significantly 


overpredicted the inside thermal resistance and led to nonsensical results for 
h,. 
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The smooth tube/twisted tape data were in close agreement with the Hong 
and Bergles (Ref. 13] correlation. Since this correlation was used by Mazzone 
(Ref. 6], his results for outside heat transfer coefficient should be fairly 
accurate for this tube. 


For the Cu/Ni finned and Korodense tubes with the twisted tape insert, the 
Hong and Bergles correlation predicted higher inside heat transfer coefficients 
than the experimental data over the entire Reynolds number range tested. 
This could have been due to a different twist ratio (for the Cu/Ni finned tube) 
and stagnant regions (for the Korodense tube). 


Experimental data for all three tubes using the HEATEX insert yielded 
- results that were well below the correlations developed and used by Mazzone 
[Ref. 6] over the entire Reynolds number range tested. 


HEATEX inserts provided the best inside heat transfer enhancement for all 
three tube. Enhancement ratios of over a factor of seven were found. 
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VIII. RECOMMENDATIONS 
Based on the results of this study, the following recommendations are made: 


Further refinement of the correlations to better account for different coolant 
inlet temperatures should be made by optimizing the Prandtl number 
exponent. 


These correlations should be used to reprocess all refrigerant condensation 
data and should be used in future refngerant condensation experiments 
conducted at NPS. Note that since these correlations were developed by 
omitting coolant velocities of 0.2 and 0.3 m/s, then the condensation data 
should do the same. 


HEATEX< inserts provide the best enhancement of tubeside heat transfer. 
Due to the data reduction technique used, HEATEX inserts should therefore 
be used in all tubeside refrigerant condensation experiments to increase the 
accuracy of the outside heat transfer coefficient. Be 


Transitional flow correlations can also be developed for the no insert 
condition since condensation experiments typically move into this flow region. 
However, since the transition region is not generally well developed, any 
developed correlations for this region will be useful only for this particular 
apparatus. 
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APPENDIX A. FLOW METER CALIBRATION 


The flowmeter used to measure flow through the instrumented tube was 
initially calibrated at 0°C and 24°C by Mazzone [Ref. 6]. However an additional 
calibration temperature of -15°C was needed for this study. The flowmeter was 
calibrated using the same procedure as Mazzone and is summarized as follows. 

The outlet of the flowmeter was disconnected from the test section and 
connected to an empty 55 gallon drum. The drum was placed on a scale with a 
resolution of 0.5 Ibs. The coolant pump was started with the flow through a bypass. 


The flowmeter was initially opened to maximum flow and then throttled back to 





10%. The weight of fluid in the drum was recorded and the timer started. The weight 
was then recorded at the end of the calibration time. This procedure was repeated 
in 10% increments up to 100% flow. The calibration run was then repeated from the 
maximum flow back to 10% flow in 10% increments. The mass flow rate was then 
calculated in kg/s based on an average of the two readings. The data for the 
flowmeter calibrations as a function of nominal flow rate are shown in Table A.1. 
A least squares linear regression was performed for the flowmeter at each of 
the three temperatures for which the calibration runs were carried out. The fmeaa 


regression equation takes the following form: 
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m = a+bN (A.1) 


where a is the y-intercept in kg/s, b is the slope in kg/s/% and N is the flowmeter 


setting in %. The results are shown in Table A.2. 


TABLE A.1. MASS FLOW RATE (KG/S) AT 15, 0 AND 24°C 


FLOW METER A 


NOMINAL FLOW MASS FLOW RATE (kg/s) 


RATE 
0°C 




























TABLE A-2. FLOWMETER CALIBRATION REGRESSION RESULTS 

Ts 
a ee 
[vinereept | 00088 | 000018 | 0007 _| 


In order to obtain flowrates at any intermediate temperatures, a simple 






interpolation procedure was incorporated into the data reduction program 
"DRPSING." The results of the calibrations for this flowmeter are shown in Figure 


A.1. 
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Figure A.1 Flowmeter Calibration. 
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APPENDIX B. HEAT LOSS TO AMBIENT EXPERIMENT 


A series of heat loss experiments were conducted to account for the amount 
of heat lost to ambient. To accomplish this, the coolant was drained from the test 
tube and the inlet and outlet capped off to prevent any air flow through the tube. 
Heater power was set at a desired level and wall temperatures. monitored by the 
DAU. Average wall temperatures were determined every minute for the first 15 
minutes and then every four minutes until an equilibrium temperature was reached 
(after approximately 60 minutes). At that point, any additional heat added could be 


considered lost to ambient. This procedure was repeated at five nominal power 





settings, 1W, 3W, 5W, 7W, and 9W. All data were taken at an ambient temperature 
of 22.0°C. 

Figure B.1 is a plot of average wall temperature versus time for the 1W and 
7W power settings. It can be seen that an equilibrium average wall temperature was 
reached after approximately 60 minutes for these two power settings (as it was for 
the other three sae: settings.) This graph shows that if the average wall 
temperature in the experiments reached 46°C or 74°C, then the heat lost to the 
atmosphere was 1W or 7W, respectively. | 

Table B.1 shows the data for the equilibrium average wall temperature as a 
function of the actual power setting (q,,,-). [he data are shown graphically in Figure 


B.2. A least squares linear regression was performed and the resulting equation was: 
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. (Fe =44.15):¢ (B.1) 
4.39WPC 


Therefore, for a given average wall temperature, the amount of heat lost to 
ambient could be determined. Equation B.1 was incorporated into the data reduction 


program "DRPSING." 


TABLE B.1 EQUILIBRIUM AVERAGE WALL TEMPERATURE RESULTS 


NOMINAL POWER | ACTUAL POWER (W) | EQUILIBRIUMT, °C) | 
=n 0.9812 46.87 
_ oa 3.122 a 
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Figure B.1 Heat loss to ambient; Average wall temperature versus time for 1W 
and 7W. 
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HEAT LOSS TO AMBIENT 
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Figure B.2 Heat loss to ambient; Average wall temperature versus heat input. 


119 


APPENDIX C. FRICTIONAL/HEAT LEAKAGE EXPERIMENT 


As coolant flows through the thermocouple wells/mrang box and along the 
tube, a small temperature rise results from frictional dissipation. When the coolant 
inlet temperature drops below ambient, an additional temperature rise can be 
attributed to heat leakage into the coolant from ambient. Since there was no means 
available to determine these terms separately, they were considered as one heat 
leakage term (4q,..;)- 

To determine this heat leakage term, the power to the heating elements was 
secured and coolant passed through the tube being tested at a given flowmeter 
setting. Any temperature rise detected was therefore due to frictional effects or heat 
leakage from ambient. This temperature rise was determined at flowmeter settings 
of 15% and 100% for each tube and insert condition at all four nominal inlet 
temperatures. 

Data showed that flowmeter setting, tube and insert condition had no 
substantial effect on the temperature rise and that only the coolant inlet temperature 
made any difference. It can therefore be concluded that the heat leakage term was 


primarily due to heat addition from ambient. The results are shown in Table C.1. 
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TABLE C.1 HEAT LEAKAGE RESULTS 


INLET TEMPERATURE TEMPERATURE RISE 





These corrections were incorporated into the data reduction program 
"DRPSING" by directly subtracting this temperature rise from the measured outlet 


temperature. 
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APPENDIX D. POWER CALIBRATION 


Actual heating element power (and therefore heat flux) was based upon the 
applied voltage and total heating element resistance. To measure voltage and 
current, the Data Acquisition Unit made use of a voltage and current sensing unit. 
However, a power factor needed to be applied to these sensed values to determine 
actual values required for the power/heat flux calculation. 

Voltage and current meters were placed in the circuit downstream of the power 
source and prior to the first heating element as shown in Figure D.1. The voltage 
rheostat was varied over its entire range and voltage and current sensed by the 
meters and in the DAU were recorded (Table D.1). For each setting, actual SEG 


was calculated by: 


P=VI (D.1) 


and total heating element resistance by: 


Roe ewe (D.2) 
ry 


[In addition, voltage (PF,) and current (PF,) power factors were determined by: 


WZ 





PF = ae and PF, = | 4cTUAL (D.3) 


: CH29 Lcz25 
where the subscripts "CH25" and "CH29" refer to the DAU sensing unit channel 
assignments and "actual" to the meter readings. All the above readings and 
calculations were done for five rheostat settings as shown in Table D.1. 
An average heater resistance (R), current (PF,) and voltage (PF,) power factors 
were then calculated. Heater power could then be calculated by one of the following 


forms of equation: 
P = (PR OPr yi. foe (D.4) 
or 
P = (PF.V..)iR (D.5) 
The current power factor varied too much to be utilized in a simple calculation. 
Therefore, it was decided to make use of equation (D.5). The final equation to 


determine the average heating element power (q,,,, in Watts) for a given rheostat 


setting was: 


Vape = (59-66 V, 


H29 


P/49.6 (D.6) 


This equation was incorporated into the data reduction program "DRPSING." 
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TABLE D.1 POWER CALIBRATION RESULTS 
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Figure D.1 Electrical schematic of power calibration experiment. 
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APPENDIX E. TUBE AND INSERT SPECIFICATIONS 


TABLE E.1 TUBE CHARACTERISTICS/SPECIFICATIONS 


PHYSICAL Cu Cu/Ni Cu/Ni 
CHARACTERISTICS SMOOTH FINNED KORODENSE 
Inside Diameter (D,) 13.26mm 10.16mm 13.35mm 
Outside Diameter (D,) 15.88mm 14.00mm 15.85mm 


T/C Diameter (D,,.) t/c 14.36mm 12.62mm 14.53mm 


Wall Thickness 


Overall Length (L) 152.4cm 153.7cm 152.7cm 
Heated Length (L,) 140.6cm 143.0cm 142.4cm 


| 

(Lr) 

[haa 
(Ls) 


T/C Location (1) 


TA «*Y«S Rm? | O0RHOm™ GG) [0.0600 


T/C Attachment Method Silver Soh Solder Tach Weld/ 
Solder Liquid Metal 


measured trom tube inlet; measured from start of heated length; (5 
finned heated length. 
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TABLE E.2 INSERT CHARACTERISTICS/SPECIFICATIONS 


thickness, 4 
width, w 
180° pitch, H 
twist ratio, y 
material 


HEATEX SMALL | LARGE 


diameter 10.16mm 13.26mm 
material stainless steel stainless steel 
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APPENDIX F. UNCERTAINTY ANALYSIS FOR THE INSIDE 
HEAT TRANSFER COEFFICIENT 

Data file SMNIO3 (smooth tube, no insert, + 10°C/10kW/m7’),was chosen to 
conduct an uncertainty analysis of the experimental determination of the inside mean 
heat transfer coefficient, h,. In order to assess the effects of coolant velocity on the 
uncertainty in the measurement of h, this uncertainty analysis was performed at 
coolant velocities of 0.4 and 1.4 m/s (for the same data set). The calculations below 
are shown for the 0.4 m/s coolant velocity. Table F.1 summarizes the results of the 
uncertainty analysis for both 0.4 m/s and 1.4 m/s coolant velocities. 

The uncertainty analysis utilizes the procedure suggested by, Kline and 


McClintock [Ref. 24]. This procedure states that if 


ING= TRACE, Xo Cage eee x) 


then the uncertainty in R, namely, dR is given by the equation: 


2 Z 2 





eee ee ee 




















where x, is the measured variable and 6x, 1s the uncertainty in the measured variable. 
The uncertainty in the coolant temperature measurements is due to uncertainty 
in the voltage measured by the thermocouples. For the thermocouples used, this 


value was assumed to be 4 microvolts or 0.1°C. The uncertainty in AT is found from 


ey 


the uncertainty in the coolant inlet and outlet temperatures, 7, and T,,, and can be 


calculated as follows: 


T. = Tas) TUS) 


Taking the partial derivatives we obtain: 


oT oT 1 


aT(14) dT(15) 2 


The uncertainty in the measurement of T,, is thus: 











05 
2 ee 
oT oT 
oT = cian tee Pag) eer urceaer, MF 
° ~ || 3Td4) aT (15) 


Based on the thermocouple uncertainty, 5T = 0.1°C, we obtain: 
dT, = 0.071°C 


The coolant inlet temperature is calculated from a single thermocouple. The 


uncertainty in this measurement is calculated as: 


T, | 
8T(13) 
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and 








zZ 
oT . 
|e |e) 
ji 0T(13) 
or, 
8T,, = 0.1°C 
Since, 





aT. i aT, ; 
dAT) ° \aaT 


then the uncertainty in AT is calculated from: 
SAT = |(sT,)*+ (87, 


Or, 


OAT = 0.1225°C 


For a measured AT of.3.7°C, the uncertainty is therefore: 


SOT _ 0.0332 = 3.32% 


AT 
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Coolant velocity is calculated from: 

77 

© pA, 

For a coolant velocity 0.4 m/s, m was calculated to be 0.05954 kg/s. The uncertainty 
in the reading of the flowmeter is defined as the scale interpolation. This value 
corresponds to one-half the value of the smallest marked increment. In this case, the 
scale interpolation factor, 5,;,, is 0.5%, which corresponds to 0.00101 kg/s resolution. 
In addition, a time-wise jitter, 5,,, of 0.5% was noted during the runs. Therefore, the 


uncertainty in m is given by: 


For a the measured coolant mass flow rate: 


GE UPS = BH 
m 


The uncertainty in the cross-sectional area of the tubes was estimated from: 


2 








where 6D; is given as 0.1 mm based on tolerances supplied from the manufacturer. 
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Therefore, for an inside diameter of 13.26 mm, 


dA 
a = 0.0107 = 1.07% 


¢ 


Assuming that the uncertainty in the evaluation of the fluid properties is negligible, 


the uncertainty in the coolant velocity is now given by: 








or for a coolant velocity of 0.4 m/s, 


bv 
= = 0.0262 = 2.62% 


¢ 





7 ee 


~The uncertainty in the calculation of coolant Reynolds number is given by: 

















-05 
SRe sv.)° (sp,\’ 
Re. - Vv yD 
or substituting appropriate numbers we find 
dRe 
= 3050282 >. 225770 
Re. 


The uncertainty in the calculation of the heat transfer rate to the coolant is 


related to the uncertainty in the measurements of coolant mass flow rate and coolant 
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temperature rise where: 


or 


24 = 0.0409 = 4.09% 
q 


The uncertainty in the heat flux is then: 








0.5 








Here, the uncertainty in the length of the tube, 5L, is assumed to be 5 mm while 5D,, 
is assumed to be 0.1 mm (again from given manufactures tolerances). Substitution 
of the appropriate numbers yields: 


5q M 
é 


It 


= 0.0147 = 1.47% 
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Finally, since h, is given by: 


It 


pm Owe 
i” IMTD 


the uncertainty in the LMTD must be calculated. The LMTD is given by: 
ici iy) ol 
T_.-T., 
ln Coe 
| 7 


In order to calculate the uncertainty in the measurement of the LMTD, the following 





equations were used: 








OLMTD _ 1 
dAT eT 
T.-T 
dLMTD AT((T,,-T..)-(T4,-T)] 
— in| ov T)(T,,-T. 
n i. 1. 
dLMTD AT 
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dLMTD _ AT 


aT, je 
| oe 


The uncertainty in the LMTD 1s therefore given by: 





y 4 
SLMTD «| ( 9LMTD sar\" | ALMTD 57 | , | SLMTD 57 OLMID sr 
AT pw aT, aT 




















Before the uncertainty in LMTD can be calculated, the uncertainty in the 
measured average inside wall temperature must be determined. The average inside 
wall was determined from a total of twelve thermocouples, three at each of the four 
longitudinal positions. The uncertainty in a single wall temperature measurement is 
due to the uncertainty in the voltage measurement of the thermocouple (65T = 0.1°C) 
and the uncertainty due to the influence caused by the various attachment techniques 
used. Obviously, the variations seen in the wall thermocouple readings were due to 
this latter uncertainty and some account has to be made of this as this is by far the 
dominant uncertainty. If the average temperature at a given longitudinal position is 
assumed to be the required temperature, then the standard deviation of the three 
thermocouple readings at that position could be assumed to be this latter uncertainty. 


Therefore, the uncertainty at a given longitudinal position can be represented by: 


134 


7 is ey 


_ I) 
ST yong) = ap 8T 


where i = 1,2,3,4 for each of the four thermocouple longitudinal locations and 


= is the standard deviation at each of these locations mentioned above. 


Since the average inside wall temperature at each position increases along the 
length of the tube, the average inside wall temperature for the entire tube must be 
thought of as the average of the four longitudinal positions’ average temperatures 


(Tos). The average inside wall temperature is then determined by: 





oT 
ST 
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Considering position 1, the average inside wall temperature is: 


ap _ T(10)+T(11)+T(12) 
i care 


pos 


For the three thermocouple readings of 
T(10) = 25.99°C, T(11) = 25.74°C and T(12) = 26.01°C, 
the average inside wall temperature for position 1 is 
ie = 25.9E 
The standard deviation at this location is given by: 
S.D. = 5T, = (ho (T,-T In)? = 0.123°C 

Therefore, 

5T. 

— Ps! 2 47x 103 = 0.47% 

pos] 


Repeating this for the other three locations gives: 


8Taa2 = 9.497°C 
ST ox3 = 9.634°C 


ST ass = 9.145°C 
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Therefore, the uncertainty in the measurement of average inside wall temperatures 


1S: 
ST -{[(0.25)(0.123)F +[(0.25)(0.497)] +[(0.25)(0.634)[ {(0.25)(0.145)F) = 0.207 


After substituting the appropriate values for dLMTD/dAT, dLMTD/0T,, 


dILMTD/dT., and d3LMTD/dT,; the uncertainty in the measurement of LMTD is: 


SLMTD = 1.122°C 


and 


OLMTD _ 1.122 = 9.0465 = 4.65% 


LMTD 24.13 


The uncertainty in the overall heat transfer coefficient, U., is then given by: ~ 


Or 


1 


h. 
= = (0.0417) + (0.0463)"|2 


rf 


dh. 
=—— = 0.0625 = 6.25% 
dh 


i 
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dh. 
— = 0.0625 = 6.25% 
dh 


i 
Table F.1 summarizes the results of the uncertainty analysis for the smooth copper 
tube with no insert for coolant velocities of 0.4 and 1.4 m/s,. respectively at a nominal 
inlet coolant temperature of +10°C. 

The uncertainty analysis neglects any uncertainty in the physical properties of 
the coolant. For the most part, this is probably a reasonable assumption. The 
exception to this may be the viscosity due to its sensitivity to temperature. Hence, the 
uncertainty reported here is likely to be conservative. 

It is curious that the relative magnitudes of the uncertainty for the low and high 
coolant velocities are nearly the same. Two different mechanisms dominate in this 
calculated uncertainty at low and high coolant velocities. In terms of the mass flow 
rate, the time-wise jitter and scale interpolation terms are the same for both coolant 
velocities. However, the actual mass flow rate is less at the lower flow rates. Hence, 
the uncertainty in the mass flow rate calculation is more significant at the lower 
coolant velocities. At higher coolant flow rates, the uncertainty in the calculation of 
LMTD becomes the dominant term. This is because the coolant temperature rise is 
smaller at high coolant velocities but the uncertainty in the temperature 


measurements remains the same. 
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TABLE F.1 SUMMARY OF RESULTS FOR SAMPLE UNCERTAINTY ANALYSIS 

















— on27o | 10.67% , 
| 
| ; 





4.09% : 10.69% 
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0.71% 


10.72% 


16.32% 


19.53% 
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1838 
ta. 
1815S |! 
1838 
1835S 
1042 
14S 
1Q@6S 
188 
1868S 
1830 
1835 
1125S 
1132 
iss 
1140 
1145 
L158 
L1S5 
1168 
117 
118 
1215 
1220 
i PaPapy 
1238 
1235 
124 
1245 
1238 
L2ss 
1268 
1261 
1263 
1264 
1265 
1266 
1267 
1268 
1269 
127 
Lz? 
L273 
1274 
1275S 
L276 
bes? 
1283 
1284 
128S 
1286 
1287 
1299 
130 
L231 
115 
L32S 
1345 
1585 


APPENDIX G. PROGRAM LISTING 


PuRBose RESINS 
QSE 1s program collacts and processas single instrumentad tube da 


CREATEO: APRIL 1, 1992 (S. MEMORY) 

SEEP 

PRINTER IS 1 

PRINT USING °4X,°°SELECT OQPTION**° 

PRINT USING °6X,°°@ TAKING OATA OR REPROCESING PREVIOUS OATA’** 
PRINT USING “6X,°°1L PURGE FILES*°” 

PRINTER IS 7@1 

INPUT Icall 

IF Icall#@ THEN CALL Main 

IF Icallel THEN CALL Purge 

ENO 

SUB Main 

COM /Ce/ Gl7) 

COM /Fld/ Ift 

COM /Nus/ [in,Tsat,Qdol,Hnus,Kf Rho? Hfg,Muf Oo, [tube 
COM /Wall/ Doa( 4) iOraca), Kma( 4), Tact ,Oroot(4),Ld(3) 
COM /Wil2/ Oelta,Isat,Nsets,Hod,Cia(4),Alpaa(3) 

OIM Em?((17),Tp(4), T7117), Hol 4) ,Qdp(4),U0( 4), Pel 4), Xstar(4),Nux(4),Nuss(4) 
OIM Kam lO) Xam 204) Xm 30 Od) eh), Mx S( 4), Xxx6(4) 
OIM Num( 4),Nuag( 4), bdd(4),Nu(4),X%3s( 4) 

OATA 1.80,8.172,0.226,0.063,5.172 

REAQ Ciale) 

OATA @.10086091 ,29727.94369,-76734S. 8295, 788255995. 81 
OATA -9247486589,6.97688Ei1,-2.66192E13,3.94078E14 
REAQD Cle) 

OATA @.019875,0.014000,2.019859,8.8159,8.8 

OATA @.013259,20.010168,8.013350,0.01346,0.2 

OATA 386.0,42.98,42.98,21.9,8.2 

OATA @.000SS88,3 

REAQ QDoa(¢),Oial*e),Kmal(e),Oelta,Hod '! HodeH/0: 

IF Itube#@ THEN 

OATA 6.76S,37.41,68.06,98.71 


REAO Ldie) y LENGTH/OIAMETER RATIO FOR 4 T/C LOCS (SMOOTH) 
L#1.486 | Heated length of smooth tube 
ENO [F 


IF Itube#l THEN 
OATA 2.05,480.880,79.55,118.29 


READ Ld(e) ! FIN LENGTH/OIAMETER RATIO FOR 4 T/C LOCS (CU/NI FIN) 
Lol s2zs ' Heated finned langth of Cu/Ni fin tube | 
ENO IF 


IF [tubee2 THEN 
OATA 7.908,38.39,63.792,95.987 . 
READ Ldie) | LENGTH/OTAMETER RATIO FOR 4 T/C LOCS (KORQOENSE ) 
Lel.422 ! Heated Length of Korodense tube 
ENO IF 
Jset 22 
Qkaceptel 
BEEP 
INPUT “ENTER MONTH, OATE ANO TIME (™M:00:HH-MM:SS)° Otgs 
QUTPUT 789: °TO°:. Otgs 
BEEP 
Ihardel 
INPUT “WANT A HAROCOPY PRINTOUT (1L*=OEF*YES,828N)°, [hard 
BEEP 
Iinel 
IF Ihardel THEN PRINTER [S 791 
BEEP 
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1360 
P65 
1370 
P3575 
1380 
FIN 

[383 
13584 


INPUT “GIVE A NAME FOR THE NEW DATA FILE” ,Files 


CREATE BOAT FileS,20 

ASSIGN @File TO Files 

BEEP 

INPUT “ENTER TUBE CODE (@=SMOOTH CU,1=FINNEO CU/NI,2=KOROOENSE, 3=TIT 
»" [tube 

BEEP 

INPUT “@2=NO INSERT, L=TWISTEO TAPE (WIO0E), 2=HEATEX (LL), 32HEATEX (5S) 


, 4"*TWISTEO TAPE ( NARROW)" Insert 


1385 
1386 
1388 
Lacs 
1330 
1391 
L522 
Pass 
[535 
1402 
1405 
1430 
143S 
1S9S 
1620 
1690S 
1612 
1615S 
1620 
L62S 
1632 
1635 
1640 
1646 
1647 
1649 
16S 
1651 
1653 
1655S 
1669 
1661 
1665 
1666 
1667 
1668 
1670 
1671 
1672 
1673 
1674 
1675 
1676 
L677 
1678 
1679 
1680 
1681 
1682 
1683 
16384 


IF Insert#l THEN Hod=4 
IF Insert#4 THEN Hod=3 
BEEP 
INPUT “NOMINAL SUMP TEMP. O=-10 C, 190 C, 2"410 C, 32+20 C*,I temp 
BEEP 
Egrat=S4 
INPUT “ENTER EG CONCENTRATION (WT PERCENT),OEFAULT = S4%",Egrat 
BEEP 
ENTER 709:0tg$ 
OUTPUT @File;:Otgs 
QUTPUT @File:Itube,Egrat ,Od1,0d2,0d3,0d4,0dS 
PRINT USING “10X,""FILE NAME: ve" 128°: Files 
PRINT 
Tout=l a 
BEEP 
INPUT “WANT TO CREATE AN OUTPUT FILE? (1=0EF*YES,O=2NQ)", lout 
IF Tout#l THEN 
BEEP 
INPUT “ENTER A NAME FOR OUTPUT FILE°,Fout3 
CREATE BOAT Fout$,§ 
ASSIGN @Fout TO Fouts 


ENO IF : 

IF [tube=@ THEN Otc#.01436 | OITAMETER TO T/C BURIAL (SMOOTH) 

IF [tube#l THEN Otc#.0126 ' OTAMETER TO T/C BURIAL (CU/NI FIN) 
IF Itube#2 THEN Otc#.01453 ' OTAMETER TO T/C BURIAL (KOROODENSE >) 


Do=Doal [tube) 
0ol=Ooa( Itube) 
O01=#O0:ia(Itube) 
KmekKma( [tube ) 


Ax2=P[001°2/4 | Cross-sectional area 

IF [tube=2 THEN Aax«=.00014 | BASED ON WOLVERINE OESIGN OATA 
Ro=PI*Doetl : 
Ar,=PTeOieL 

IF Itube#2 THEN Ao=.0SeL | BASEO ON WOLVERINE OESIGN OATA 


IF Itube=2 THEN Ar*A0/1.186 |! BASEO ON WOLVERINE OESIGN OATA 
Per im2=PI «01 

Anarea®Ple(Q0°2-01°2 7/4, 

Ldd( l)=L/(4, 01) 

ads yal (2 seo) 

Ldd(3)83.°eL/( 4.01 ) 

Ldd(4)eL/O1 

IF [tube#0 THEN 


Ulent=.Q627 | UNHEATED INLET ENO OF COPPER SMOOTH TUBE ONLY 
Ulen2*.0S72 | UNHEATEO DUTLET ENO OF COPPER S001 ra) eee ona: 
ENO IF 

IF [tube*l THEN 

UlLenl=.045 | UNHEATEO INLET ENO OF CU/NI FIN TUBE ONLY 
Ulen2=.06 UNHEATEO OUTLET ENO OF CUYNI EIN: TUBE  GNEY 


Ocuo#.01587S |! UNHEATEO OUTSIOE OIAMETER OF CU/NI FIN TUBE 
Ocupet.015259 !' UNHEATEQ INSIOE CDAMETER CP eCUy Nie er 
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1685 
1686 
1687 
1688 
16893 
1692 
i631 
Hee 
1695 
16934 
1700 ! 
17@1 
1782 
17@3 


/3)?”" Vegst 


1704 
1705 
1706 
1707 
1708 
1709 
1710 
Lele t 
ie i 2 
R713 
1714 
i715 
1716 
ei? 
1718 
iis 
1720 
el 
Net 2 
23 
1724 
eS 
26 
Vr eeage 
1728 
ied 
1730 
Uy 
733 


beugevtl2 
Leuo#.097 


Perim=PI eQcu1 
Anarea=Pl«(Qcu0*l-Ocu1*2 )/4. 


ENO IF 
IF I[tubes2 THEN 


Ulenl=.050 
Ulen2=.055 


ENO IF 


Hr =5Q 


RRINTER oe! 


BEEP 
INPUT 


SMe A EOmINEET SMOOTH END ORV EU/NT FIN TUBE 


HEATEQ OUTLET SMOOTH ENO OF CU/NI FIN 


UNHEATEQ INLET ENQ OF CU/NI KOROQENSE TUBE ONLY 
UNHEATEO OUTLET ENO OF CU/NI KQROOENSE TUBE ONLY 


| HEATER RESISTANCE 


“WHAT IS REQUIREQ NOMINAL COOLANT VELOCITY IN m/s (MAX. 


OUTPUT 709: "AR AF1L3 AL1L3 VRS” 
OUTPUT 709; "AS SA“ 
Vsum=@ 
FOR J#1 TO S 

ENTER 709;E€ 

VsumeVsumee 
NEXT J 
Emf(13)98Vsum/5 
TCL3 *FNTvsv( Em((13)) 

Tavgl=T(13) 

IF Tavgl<-15 OR Tavgl>24 THEN 


PRINT “COOLANT TEMP. NOT IN RANGE -15 TO 24 C" 


GOTO 1703 


ENO IF 


Rhoegl =FNRhoeg( Tavg! ,Egrat ) 
Mdotl=Rhoeg! *Ax*Vegst 
PealS=(Mdotl+.0098 )/.0019 

PceaQ=(Mdotl+.00414)/.00205 
PeaZ4=(Mdotl+.00329)/.00225 
IF Tavgl<0@ THEN 


Po(Q)=PcalS+( ( Tavgl +15 )9(Pead—-PealS) )/15 


BESe 


Po(Q)=PcaQ+( Tavgl*(Pea24-Pcad))/24 


ENO IF 


PRINT “SET FLOWMETER REAQINGS CORRESPONOING TO: “ 


PRINT USING “6X,°°% OF METER A = “",Q00.0";:Pc(O) 


PRINT 


PRINT “HIT CONTINUE WHEN READY” - 


PAUSE 


1734 Repeat: 
PRINTER [IS 7@1 


1735 | 
1737 
1738 
"1739 
1742 
1741 
1742 
1744 
1745 | 
1746 
1747 
1748 | 
1749 
1750 | 
17S1 


Ido#l 


ONGKEY -O.15 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
Ido = 
BEEP 

INPUT 


USING 
USING 
USING 
USING 
USING 


RECOVER 1734 
Sox, See | 
"6X," "@=TAKE 


OPTION = 
CA Ae 


TUBE 


"6X,""L#SET Heater Tape Heat Flux (QEFAULT ) 


“GX, ""2=CHECK FLOWMETER % SETTING 


"4X," "NOTE: 


desired ootion 


Ido 


Set default value for in 
IF Ido>2 THEN [do#2 

Take data option 

IF Idoz@ THEN 1816 


KEY O =2ESCArE™ 


put 
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m- 4 


= 1.4m 


IF {do=#=2 THEN 17@3 
Loop to set heat flux and check temperatures 
IF Ido#l THEN 

INPUT “ENTER DESIRED HEAT FLUX" ,Qdpwan 


PRINT USING “4X,"°"* Qdp (wanted) Qdo (actual) Tin Geli 


PRINT USING “4X,"%"  (W/m*2) (W/m*2) (C) (C) 
Read thermocouple voltages 
DUTPUT 709; "AR AFL ALLS VRS" 


FOR [#1 


Lots 


OUTPUT 709: "AS SA” 
Vsum29 
FOR J#l TO 2@ 


ENTER 7@9;E 
Vsum=VsumeE 


NEXT J 
Emf( I )#Vsum/2Q 
TCL) *#FNTvsv(Emf¢ I )) 


NEXT I 


Compute average temperature of coolant outlet 
Tout#(T(14)4T(15))/2 
Delt=Tout-T( 13) 
OUTPUT 709; "AR AF2S ALZS VRS" 
OUTPUT 7@9;"°AS SA” 
ENTER 709; Amp 
OUTPUT 709: "AR AF2Z9 ALZY VRS” 
OUTPUT 709: "AS SA“ 
ENTER. 7@9;Volt 
Qtape=(SQsVolt )°2/Hr 
Qdptape=Qtape/Ao 


PRINT 


PRINT USING °4X,S(MZ.30E,3X)":Qdpwan, Qdptape,T(15),Delt,Qtape 


PRINT 


PRINT USING 


PRINT 


PRINT USING 


PRINT 


PRINT USING “4X,° "TEMPS Al POSITION 3° >. S(QO0700 (5x )°7 1 4 aes eee 


PRINT 


PRINT USING 


WAIT 2 


GOTO 1762 


ENG TF 


TAKE OATA IF I[m2Q@ LOOP 


Twm2Q 
Twi2@ 
PRINTER [IS 


OUTPUT 7@9:; 


FOR [1 TO 
OUTPUT 
Vsum=Q 


701 

"AR Ab ARTS VRS” 
1s 

US ees icl=) S015 


FOR J=1 TO 2@ 
ENTER 7@09;E 
Vsum2Vsum+E 


NEXT J 


Emf( I )#VUsum/2@ 


NEXT I 
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“ 


ab 


"4X, °° TEMPS AT POSITION 1°” ,SCOOU 200. 3X) ; 1 1G) ea erie 


"8X, °°TEMPS AT POSITION 2°" , 3°00. 00.35X)": 1° 7 ice ee 


"4X, “"TEMPS AT POSITION 4°*,3(000.00,3X)";T(1),T(2),T(3 


DATA ANALYS 


FOR T= alo 


is 


1S 


TCI )*FNTvsv(Emf¢( I) ) 


NEXT I 

FOR [#1 TO 
TwitTwit 

NEXT I 

TwiaveTwi/l 


te 


1a 


ees 


Touts(T(14)4T(15))/2 


Delt=Tout-T 
Hlc=Q. 

IF I[tube=2 

IF I[temp=2 

IF Itemp=l 

IF Itemp=2 

IF [temp=3 

DeltsDelt-H 
ENO [F 

Tout *Tout-H 


OUTPUT 709; °AS SA“ 


(13) 


| Heat leakage correction 


THEN 


THEN Hic#.18 
THEN Hlce#=.Q8 
THEN Hlc#.07 


THEN Hlc=2. 


le 


Le 


ENTER 729: Amp 
OUTPUT 709: °AR AFZ29 ALZ29 VRS" 


OUTPUT 709;°AS SA* 


ENTER 709:Volt 
Qtapes(GO*Volt )*2/Hr 


Qdptape=Qta 
Jsetetset+l 


pe/Ao 


IF Okaccpt#=@ THEN Jset#Jset-1l 


PRINT 

PRINT USING 
PRINT 

IF [tube=@ 

IF I[tube#2 

IF I[tube=l 

IF [temp 

IF [temp=1 

IF [temp=2 

IF Itemp=3 

IF Insert= 
IF Insert=l 
IF Insert=4 
IF Insert#e 
IF Insert=3 
PRINT USING 
PRINT USING 
PRINT 


TavgsT(l3)+0elt*.S 


“10X,""Data set 


THEN 
THEN 
THEN 
THEN 
THEN 
THEN 
THEN 


THEN, PRINT USING 
PRINT USING 
THEN PRINT USING 
PRINT USING 
THEN PRINT USING 


THEN 


THEN 


"LOX 


“10X," "Nominal set coolant velocity 


PRINT 
PRINT 
PRINT 
PRINT 


“PRINT 


PRINT 
PRINT 


number 


HEAT LEAKAGE CORRECTION AT -1@ C 
HEAT LEAKAGE CORRECTION AT @ C 

HEAT LEAKAGE CORRECTION AT +190 C 
NO HEAT LEAKAGE CORRECTION AT +20 C 


Pee i Set 


USING “10X,""°Smooth tube”"”” 

USING “1LOX,°"*Korodense tube""" 

USING “L@xX,° “Cepper/Nickel fim ~~ 

USING “10X," "Nominal sump temp. = -10 C*"" 
USING “L@X,""Nominal sump temp. = OC 
USING “10X,"“"“Nominal sump temp. = +190 C°*" 
USING “L@X,""Nominal sump temp. = +20 C°"" 


~” "Nominal set meat flux ~# 


Rhoeg*FNRhoeg( Tavg,Egrat ) 
Nuueg#FNNueg(Tavg,Egrat ) 
Mueg*Nuueg *Rhoeg 
Coeg*FNCpeg( Tavg,Egrat ) 
Keg#FNKeg( Tavg,Egrat ) 


Preg=Cpeg*Mueg/Keg 


Mdot =FNFmcal(T(13),Pe(Q@)) 
Veg=Mdot/( Rhaeg*A» ) 
Reeg=Veg*01/Nuueg 


Res=4eMdot/(PI *Muege(01-4*0elta)) 


Qdot=Mdot*Cpoeg*Delt 
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"“10X,"°"No insert used"”"” 
“10X,*"Wide twisted tape 
“LOX,""Narrow twisted tape insert used" ”” 
"10X,“"“Large HEATEX 
"10X,“"Small HEATEX insert used"”” 


Lnsert used~”“ 
insert used"”" 


"",%2Z.30E" :Qaptape 
-2eeO. 0 | Vegst 


2048! [ITERATE TO FINO SINS TOE WAEE En eRa URES 


2049 Tml=Twiav 

2252 Tm2=Twiav 

2@S1 Tmla( TmitTm2)/2. 

28S2 Ctr=Qdot*L0G( 0te/01)/(2. Pl eLekKm) 

28S3 Tm2*2Twiav-Ct1 

2@S6 IF ABS( Tm2-Tml)>.@@1 THEN GOTO 2@Sl 

20S7 FOR [=l1 TO l2 

20S8 TCL )#T¢( lL )-€ti 

2839 Tuwm=TwmtT([) 

2060 NEXT [ 

2061 TwmeTum/ 12 

2062 Nuuegw®F NNueg( Tum, Egrat ) 

2063 Muegu=Nuuegu*Rhoeg 

2064 Twml ac T O12 +TO LL +7618) )/3 

c@6S Twm2=(7T(69)4+T(8)+T(7) 0/3 

<@66 Twm3=( 766 )47T(5)4+T( 4) )/3 

2067 TwmS2(7T(3)4T62 4TC 1) 9/3 

2268 Twmi2acTCL2 TCL +T¢ 18 )+T69)4+TC8 +767) )/6 

2269 Twm loser TOL2 +TOL lL +70 100470904718 4T07 4T1G 4T(S 471 4))/9 
22972 Tim=( Tout-T(13))/LOG¢( Twm-T( 135) )/¢ Twm-Tout )) 

2071 Qloss=( Tum-44.15)/4.39 |! HEAT LOSS TO ATMOSPHERE 

2272 IF Twm<4S THEN Qloss=2 

2073 Qtape=Qtape-Qloss 

2074 Qdpr.=Qdot/A1 

2@7S AlfaizQdp1/Tlm 

2276 |! FIN CORRECTION TO ACCOUNT FOR CONOUCTION TO ENOS OF TUBE 
277 FPL a(Km*Perim*Anarea)*.S 

2978 Alfairl#Alfan i Css 
20793 AlfarcsAlfail 

2282 Alamla( (Al far2*Perim)/(Km*Anarea))*.S 

2@81 FF3eFFLeALPai2* .Se( Tum-T(13)) 

2282 FF4eF Fl eAlfar2*.S*( Tum-Tout ) 

2@83 Func6 =F f3*FNTanh(Alaml*Ulentl ) 

22084 Func7#F f4eFNTanh(Alaml*Ulenc ) 

2@8S Func82Al fai2Z*A1*( Tum-Tavg) 

2@86 Funca*Func&+Func7+Func8-Qdot 

2287 Ofunc62( .S*FuncG/Alfaiz del. oF F3ZeUlenl/( 1. +FNCosn(l.*Alaml*Ulen] )) 
2988 Ofunc7#( .S*Func?/Al fairc +2. °F F4eULen2/(1.+FNCosh( 2. *Alaml*Ulenc )) 
-@83 Ofunc8#A1*( Tum-Tavg) 

2-2@9@ Ofuncx#O0func&+Ofunc?+Ofuncs 

20a Alfarl#Al farl-Fune«}/Ofunex 

2092 IF ABS(Al farl-Alfai2)>.@S THEN GOTO 2079 

2ege * AlfarewAlfarl 

2045 Cfeg#(Mueg/Muegw)*.14 

2096 | FINO TEMPERATURE INCREASE OUE TO UNHEATEO LENGTHS 

2297 Oeltin=Func6/‘Mdot *Caeg) 

2938 ° Deltout=Func7?/(Mdot*Cpeg) 

Zod Trlin=T(13)+0eltin 

2122 Thlout#Tout-Oeltout 

21g Delhl=( Thlout-Thlin) 

21@S IF [tube#l THEN 

21@6 Heatin=eicul*Oelhl/(L+LecurtlLcuo) !' HEATEO SMOOTH INLET SECTIGN FOR CU 
/NI FIN TUBE 

2127 Heatout=LcuoeDelhl/(LtLcurtlLcuo) |! HEATEO SMOOTH OUTLET SECTION FOR 
CU/NI FIN TUBE 

21@4 ThlineTrlintHeatin 

etl Thlout#Thlout-Heatout 

a Delhl=(Thlout-Thlin) 

zis ENO IF 


2415 Qhl=Mdot *Cpeg*Delhl 


eis Qdphi=Qni/Ar 
Bali? Qdp«x=Func8/A1 
21198 Teol#Thlint(Oelhl/4) 
e119 Teo2=Thlint(Oelhi/2) 
2129 Tco3#Thlint( 3eDelhl/4) 
Zl2i TtelsT(t3)+(Ld(@)*O1*Delrl/L) 
Bec Tte2=T( 135 )+( idl teOr1*DelhAl/L) 
2123 Tee3eT(13)4+(bd(2)901*Oelhl/L) 
2124 Tte4aT(13)4(bd(3)eO1*eDelrl/L?) 
2126 Tlmall#( Thlout-Thlin)/LOG( ( Tum- Ena ee eee Alt (1.2 sees) 
Zle7 Tlmal#(Tcol-Thlin)/LOG( ( Tuwmt-Thlin)/( Tumi-Teol )) 1 SECTION 1 
etec Tlm22( Tco2-Tcol )/LOG( ( Tum2-Teol )/( Tum2-Tco2)) | SECTION 2 
2129 Tlm3#( Tceo3-Tco2 )/LOG( ( Tum3-Teo2 )/( Tum3-Teo3) ) ' SECTION 3 
2150 Tlm4=( Thlout-Teo3 )/LOG( ( Tum4-Teo3 )/( Twm4-Thlout )) 1 SECTION 4 
gol Tlmi2=( Teo2-Thlin)/LOG(( Tuml2-Thlin)/(Twml2-Teo? ) ) SECTION faz 
else Tlmi23=( Tco3-Thlin)/LOG( ( Tumli3s-Thlin)/(Twml23-Teo3))! SECTION 1,2,2 
Z1L35 Alpall#Qdpohl/Tlmall 
2134 Alpl*Qdohl/Tlmi 
2135S Alo2=Qdohl/Tlm2 
2136 Alp3*Qdahl/T1lm3 
Zus7 Alp4=Qdphl/T1lm4 
2138 Allocit#Qdphl/( Tumi-Ttel) 
“BSS Alloc2*Qdphl/( Tum2-Tte2) 
2140 Alloc3=#Qdphl/( Tum3-Ttc3) 
2141 Alloc4#Qdpahl/( Tum4-Ttc4) 
2143 Aloit2=Qdphl/Tlml2 
ae Te Alpti3“Qdphl/Tlmt23 
2145 Kegl#FNKeg((TcoltThlin)/2,Egrat ) - 
2146 Kego=FNKeg(( Tco2+Tcol)/2,Egrat) 
2147 Keg3sFNK eg( (Tco3t+Tco2 )/2,Egrat ) 
2148 Keg4=FNKeg( ( Thlout+Tco3)/2,Egrat ) 
2149 Kegl2=FNKeg( ( Tco2+Thlin)/2,Egrat ) 
2150 Keg lo3sFNKeg( (Tco3+Thlin)/2,Egrat ) 
Bol Numall#Alpall*O01/Keg 
2152 Nul#Alp!*O1/Keg! 
2153 Nu2*#Alp2*01/Keg2 
2154 Nu3#Alp3¢01i/Kegl 
A USS Nu4#Alp4e01/Keg4 
2 US| S Nulocl#Alloci*O1/Kegl 
2is7 _ NulocZsAlloc2*O1/Keg2 
21$8 Nuloc3#Alloc3e01/Keg3 
z2is9g Nuloc4#Alloc4*01/Keg4 
2161 Nul2*Alsl2*O1/Keg!2 
2162 NulcsS=Alpl2seO1/Kegli23 
2163 |! LOCAL NUSSELT FROM SHAH/LONOON PAGE [28, EQN 240-240 
2164 |! CALCULATE XSTAR FIRST FOR EACH THERMOCOUPLE LOCATION 
cl6S ! OEFINE COOLANT PECLET NUMBER 
-166 Peag=Preg*Reeg 
2167 Xs(t)eidd(t)/Peeg |! FOR MEAN VALUES OF A 
2168 Xs(2)=Ldd(2)/Peeg 
2169 Xs(3)=Ldd(3)/Peeg 
2172 Xs(4)sLdd(4)/Peeg 
217i FOR [=!1 TO 4 . 
Dl We. Xstar([)aid([-1)/Peeg !'§ FOR LOCAL VALUES AT T/C POSITIONS 
2175 Nuxw CP =t.502*xXstar< fy (-1/3)—.5 
Ao. | DEVELOPING HYORO ANO THERMAL BOUNOARY LAYER (S/L P148, EQN. 253) 
2176 Reamer eC OLO/P Teka tars 1 ))°C=10/9) ba0o75) 
2177 Kala(L+(Preg/.0207)*(2/3))7.5 
2178 Xx ZaPlI/CLi1S.2eXstar(1)) 
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rae Kade lL $OXkS/ (Xx 29h x1) C575) ea 


218@ XxSaC 14¢((22Q0/PI)*Xstar(1I))*(-10/9))*(3/18) 

216! Nuss(I )#(X«405,3640X«5 )-] 

2188 IF Xstar(I)<.@@0Q@S THEN 

2189 Nux(I)#l1,302¢Xstar(I)*(-1/3)-1 

219 ENO IF 

2194 IF Xstar(I)>.@@15S THEN 

2045 Nu«( I )2#4.564+8.68( 10Q0*Xstar(1I))*(-.5Q06 )*EXP(-41*Xstar(I)) 

2196 ENO IF 

22202 NEXT I 

2202 IF Insert=@ THEN 

2203 PRINT USING “3X,°" LOCAL NUSSELT VALUES"°*”" 

2.0S PRINT USING °3XxX,°° Exp. Nux(TLM) Exp. Nux(8ULK ) Pred, Nua(l ) 
Pred. Nux(2) Xstar(tc)*”*” 

2206 PRINT 

2207 PRINT USING °3X,°°POS 1 ©" ,S(MZ.30E,5X)°> Nut Nulocel,Nux( t),Nuss (ei 
star(1) 

2208 PRINT USING °3X,°°POS 2 °° ,S(MZ.30E,SX )°:Nu2,Nuloce ,Nux(2),Nuss(2),X 
star(2) 

2209 PRINT USING °3X,°°POS 3 °",S(MZ.30E,SX )°i:Nu3,Nuloc3 ,Nux(3),Nuss(3),X 
star(3) 

2210 PRINT USING “3X,°°POS 4 °° ,S(MZ.30E,SX )°;Nu4,Nuloc4,Nux(4),Nuss(4),X 
star(4) 

Zeit PRINT 

Ze ELSE 

2216 PRINT USING °3X,°°POSITION &% ‘= Exp. Nux (TLM) Exo. Nux (8ULK) 

Xstar at tc pos.””” 

Lela PRINT 

Zele PRINT USING °3X,°**POSITION 1 (LOCAL) = °",3(MZ.30E,5xX )":Nul,Nuloel,XxX 
star (l) 

2213 PRINT USING °3X,°”"°POSITION 2 (LOCAL) = °° ,3(MZ.30E, SX )*:Nu2,Nuloce,Xx 
star(2) 

cece PRINT USING °3X,°"POSITION 3 (LOCAL) = °*,3(MZ.30E,SX)"*:Nu3,Nuloc3,X 
star(3) 

Farce | PRINT USING °3X,°"POSITION 4 (LOCAL) = °",3(MZ.30E,5X )"iNu4S,Nuloc4,Xx 
star(4) 

Ree ENO IF 

Fa Pig PRINT 

wee FOR I=1 TO 4 

Zeon IF Insert=@ THEN 

Fag | Bol=.Q@668*Reeg*Preg/Ldd( I) 

2254 Bo221+.Q4¢( ReegePreg/Ldd(I))*.6666. 

Leos Nueg( I )#3.66+(861/Bb2) 

2eob ENO IF 

2250 IF Insert=l OR Insert=4 THEN 

Zoos Nueg( [)2Cira( Insert )e(1+S.484E-3ePreg”.7*(Res/Hod)*1.25)°.5 

2240 -Nuoth=™.383*Preg”.35*(Res/Hod)*.622 

cet ENO IF 

cade IF Insert=2 THEN 

eos Nueg( I j#Cia( Insert )*(Reeg”.65)*Preg”. 46 

2244 ENO IF 

2245S IF Insert=#3 THEN 

2246 Nueg( I )*#Cra( Insert )*(Reeg*.76)*Preg*. 46 

2247 ENO IF ° 

2248 BEEP 

249! CALCULATED VALUE OF INSIOE HT. TRANSFER COEFF. FROM CORRELATIONS 

2250 Hi=Nueg(4)*eKeg/01 

Zeal MEAN NUSSELT NUMBER FROM SHAH/LONOON PAGE 128, EQN 2459-246 

cone IF Xs(1I)%.03 THEN 

e255 Num( I )#1,9S3eXs(I)*(-1/3) 
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e250 ELSE 


Z2L58 Num( I )=4.3564+( .@722/Xs(I)) 

2259 ENO IF 

2262 NEXT I 

2263 PRINT USING “3xX,"“"POSITION & = Exp. mean Nu Pred. mean Nu 

Xstar at section"”” 

2264 PRINT 

2286 PRINT USING “3X,""POSITION | (MEAN) = "" 3CMZ.30E,S4)":Nul Num Ll) 

Xs(h) 

2267 PRINT USING “3X,°"POSITION 12 (MEAN) = "" 3(MZ.30€ .SX%)"° i: Nul2,Num( 2 

»,Xs(2) 

2268 PRINT USING “3X,""POSITION 123 (MEAN) = °° ,3¢MZ.30E,5X)*;Nul23,Num( 

S).X%3(3) 

2269 PRINT USING “3X,""POSITION 1234 (MEAN) = "°,3(MZ.30E,S%)"sNumall,Num 

(4),Xs(4) 

2272 PRINT 

Zeal PRINT 

Bere PRINT USING “10X,°°Mass flow rate (kg/s) eee enl. SUE i Mdot 

Ze7s! PRINT USING “1@X,""Inside Tube Dia. (m) 2 ""MZ.30E" Oral Itube) 

2274 IF Insert=l THEN PRINT USING °10X,°"18@ OEG OVER Oia. (HOO) = "", 

fe. o0e” :Hod 

2273S PRINT USING “1OX,""Inlet temperature (C) Seon e. Soe . fle, 

Be76 IF Insert#l OR Insert=4 THEN 

2277 PRINT USING “10X,""ODELTA Tape Thickness eee ic, SOE. : ela 

2278 ENO IF 

2279 PRINT USING °10X,°"OELT temp Oif. (C) =e te. oUe 7Oelt 

2280 PRINT USING “10X,°°H. T. to coolant (W) “eee. SOE + Qdot 

aol PRINT USING “10%, "cHe= |. eomeneater tape (W) = ~~ ,MZ.30E": Qtape 

Ezgc PRINT USING “1Q@X,""Prandt! number = “* .MZ.30€";Preg a ee 
== 7742 es PRINT USING “L@X,""“Reynolds number = "" MZ.30E";:Reeg 

2284 IF Insert#l OR Insert=4 THEN-PRINT USING “1X,” "Reynolds number H 

8B S = "" MZ.30E";: Res 

2235 PRINT USING “10X,""Hi from correl. (W/m°2.K) = "" MZ.50E".H1 

2286 PRINT USING “10X,""Hi from heat bal. (W/m*2.K) = “" ,MZ.S0E";Alfanr 

2 

Zeal IF Insert#@ THEN PRINT USING “10X,""Nu for no insert (Hausen) 

".MZ.30DE" ;Nueg( 4) 

2288 IF Insert#!1 OR’ Insert=4 THEN PRINT USING “L@X,”""Nu for twisted ta 

pe (Bergles2)"",MZ.350E" i Nuoth 

acd IF Insert=l1 OR Insert=4 THEN PRINT USING "1@X,""Nu for twisted ta 

pe (Bergles3)°",MZ.30E" i Nueg( 4) 

2290 IF Insert#2 THEN PRINT USING "1@X,""Nu for large Heatex (Mazzone) 

"" MZ.30E" :Nueg( 4) 

2291 IF Insert=3 THEN PRINT USING “1@X,""Nu for small Heates (Mazzone) 

"" MZ.30E”" :Nueg( 4) 

ez. IF Insert=@ THEN PRINT USING "“1LOX,""Nu from Shah/London ae 

MZ.30E" :Num( 4) 

e295 PRINT USING “10X,""Xstar for heated length = "",M2.2D0E".X%s(4). 

ee 4 PRINT USING "L@A,""E-perimental Nusselt no.= "',M2.I50E" .Numal! 

e245 PRINT 

2238 PRimUsSINGuws4, Watt t/o L2,Ll,1@° (oes... 1) (CC) = °",35(098.00,5Xx 

Relay rcniy. TC 1e) 

Bag! PRINT USING “3X,""Wall t/e 9,8,7 Geos weGo) = "(00000 ,.o% 

ie 9) eB), 147) 

ee Sys PRINT USING “3X,""Wall t/e 6,5,4 Gposseommcc) = °*,5(0D0200,2x 

me CS) CS). 104) 

2299 PRINT USING “3X, “Wall t/e 3,2,1 (gose¢)-(C) = °”, §°0@08.90,5x 

mee 3 eee), TC 1) 

2302 PRINT USING "3X,""Av. inside wall temp (C) = "",MZ.S0E". Tum 

2301 PRINT 
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2302 Nu( lt )#Nul 


2303 Nu(2)*@Nul2 
2304 Nu(3)*Nul23 
2305 Nu(4)sNumall 
2307 IF Insert=@ THEN 
2308 FOR [=i TO 4 
2312 Xxal( IT )sLOG(Nu( I) ) 
23h) Xa242(T )*LOG(Xs(I)) 
coke X<«x3(1T)*LOG(Nueg( I)) 
Zol3 Xxx 4(I )sLOG(Num( 1 )) 
2314 NEXT I 
2al15 PRINT USING °3X,°”"°MEAN VALUES AT POS 1 POS 1 
2 POSmrce POS 1234°"" 
2ol7 PRINT USING "&X,°°EXP. LOG(Nu) = °° 4(MZ.30E, 4X)" i: XxxlC lox 
xx162),Xxx103),Xxcal(d) 
Zoe PRINT USING "6X,°°COR. LOG(Nu) (H) = *" 4(MZ. 30E , 4X 1° 7 Kee 
xxX3(2),XA43(3),XxK5(4) 
Za. PRINT USING "6&X,"°COR. LOG(Nu) (S/L) = °* ,4(MZ.30E,4% "i Xxn401),X 
xx4(2),Xxx4(3),Xxx4(4) 
2320 PRINT USING °6X,°"°LOG(Xstar) oo" 4(MZ. SOE, 4X)": Xxx2 cee 
xx2(62),XKK2(3),Xxx2(4) 
ese ENO IF 
Z5e0 IF Insert#2 OR Insert=3 THEN 
2324 FOR I=l TO 4 
wo29 Xx«l( I )*LOG(Nu( I )/Preg*. 46) 
2326 Xxx2( I )slOG(Reeg ) 
coce Xx«43(1 )SLOG(Nueg( I )/Preg*.46) 
2328 NEX{ lo ee 
Zoes PRINT USING °3X,""MEAN VALUES AT POS 1 POS e 
2 POS 123 POS le34 28 
2351 PRINT USING °6X," "EXP. LOG(Nu/Pr*.46) = *" ,d(M2.350E, 4X)": Xxxl(l), 
Rex lCZ Xxx las kx) C49 
AB SIR PRINT USING "&X,""COR. LOG(Nu/Pr*.48) = “"4(MZ.30E,4X "i Xax5(l), 
Xx430 2), XK45035) , X245048) 
2334 PRINT USING “6X,°"LOG(Re) = **  4(MZ. SOE 4X “i Xx Glo 
Kee ), Xue eG), kxxe 3s) 
2555 ENO IF 
25S IF Insert=l OR Insert=4 THEN 
goo 7 FOR [=1 TO 4 
2538 X.«1¢T)*#LOG(Nu( I )/Preg*.sS) 
aoo9 Xx«2([)sLOG(Res/Hod) 
=~ 340 Xx43(1)*LOG(Nuoth/Preg*.35) 
2341 Xxx 4( 1 )=LOG(Nu( I )) 
Saae a: X45( [ )#@lOG( Preg*( Res/Hod)*1.78) 
2343 Xxx6( 1 )*#LOG(Nueg( I)) 
2344 NEAT I 
2345 PRINT USING *3X,°*"MEAN VALUES AT POS 1 PQS 1 
Ri POS shoe POS lec4d. 
2346 PRINT USING “6X,° "EXP. LOG(Nu/Pe° 155)" ene, SOE. 44 eee ee 
XM ALCS) XRG) Mxan la 
2347 PRINT USING "6X,""°COR. LOG(Nu/Pr*.35S5) = “",40MZ.350E, 4X)", Xax Sl), 
RAn5CE 1, X4KS 03 hag 4) 
2348 PRINT USING “6X,"*“LOG(Res/y) a m 40MZ.3S0E 4%)" shen oe 
Meals WX xne5 keel Cd) . 
Zoa3 PRINT 
22350 PRINT USING “6X,""EXP. LOG(Nu) = mC Cg De bake eked Cale 
Aw x42), ker8(5) koe d( 4) 
2251 PRINT USING “6X,""COR. LOG(Nu) = a TSCM Z2 SOE aAX Ioan Ole 
Xx46(2),%446(3),Xa26(4) 
Zo5e PRINT USING “6X,"“LOG( Pr (Res/y)*1.78) = “*,4(M2.30E, 4X)". XanS( 1), 
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Maes e ),Xx*%5(3),X««5(4) 


2355 


ENO IF 

IF Insert=@ ANO Reeg?2500 THEN 

Nutu#.Q@27*Reeg” .8*Preg* 3333 *Cfeg 

PRINT USING “6X,""TURBULENT S/T CORRELATION = “*,4(MZ.350E,4X)";Nu 


Xxx 7sLOG(Nutu ) 
PRINT USING “6X,""TURB. COR. LOG(Nu) ee 6 44MNZ SOE, 4X)” i Xxne 
ENO IF : 
Okaccpt=l 
BEEP 
INPUT "OK TO ACCEPT THIS SET (1L=0EFAULT=YES, @=NO)?" ,Qkacept 
IF Okacept=l THEN QUTPUT OFile;:Pe(*),Emf(*),To) 
IF Im=Q@ THEN 
Okrpt=l 
8EEP 
INPUT “WILL THERE 8& ANOTHER RUN (12YES20EFAULT ,Q2=NQ)" ,Okrpt 
IF Okrpt=l THEN 17@1 
ELSE 
IF Jset<Nsets THEN 17921 
ENO IF 
ASSIGN @File TO #- 
IF ITout#l THEN ASSIGN @File TO » 
SUBENO 
OEF FNGrad(T) 
Grad=-3:8778S7E-S-294. 71428S7E-8 eT 
RETURN Grad 
FNENQO 
OEF FNKeu(T) 
OFHC COPPER 25@ TO 300 K 
TkeT+#273.15 
K=434-,112°Tk 
RETURN K 
FNENO 
OEF FNNueg(Tec,Egr) 
RANGE OF VALIOITY: <-20@ TO 20 OEG C 
TksTet+275.15 
Nul=7. L196SQ7E-3-Tk#( 7.486334 7E-S-Tk #(2.6294943E-7-Tk*3.08335529E-10)) 
Nul24.9237638E-3-The( 4. S21L359L2ZE-S-Tk o( 1. 64357S34E-7-Thel . 3352535237 1E-12) ) 
Nu328.6S86293E-3-Tke(8.88379Q0LE-S-Tk (35. Q049S052E-7-Th #3.4999996E-10)) 
AZ=(NuZ—2 *Nuc+Nut 1/200 
AL=(Nu2-Nul~94@eA2 )/12 
AQ@*Nul-42*°AL~-17649A2. 
NuzAQ+tEgre( AL +Egr Az ) 
RETURN Nu 
FNENO 
DEF FNCopeg(Tc,Egr) 
RANGE OF VALIDITY @ TO 20 DEG C 
Tk2Te*+2735.1§ 
Col=1.G6701SSQ0E+3+Tk 96.3 
Co221.47481°SE+5+Tk 96.25 
Cp5=9.S8Q00SQQ0E+2+Tk*e7.5 
AZ=(Co3-2°Coc+Col )/200 
AL=(Coc-Cpl-900*AZ )/19 
AQ@=Co1l-40°AL-LEQQ*AZ 
CozAQ+Egre(AL+Egr Ac ) 
RETURN Co 
FNENO 
DEF FNRhoeg(T,Egr ) 
Rol=1.Q06070935E+3-T*(3.70312835E-1+T*4.0837183E-3) 


Lol 


Ro2=1.0748272E+3-Te(4.426619SE-1+T*4.0939706E-3 ) 
Ro3=1.0885934E+2-Te(S5. 735S6S5E-1+T*6.1281405E-3 ) 
AZ=(Ro3-2¢Rol+RoOl 1/200 
Al=(Ro2-Rol-980*AS 1/12 


AQ=ROl-40"eAL-1600eAZ 
Ro=AQ+Egre(AL+EgreAc ) 
RETURN Ro 

FNENO 

DEF FNPreg(T,Egr) 


Pr=FNCpeg(T ,Egr )*FNNueg( T,Egr )*FNRoeg(T,Egr)/FNKeg(T,Egr) 


RETURN Pr 

FNEND 

DEF FNKeg(Tc,Egr ) 
RANGE OF VALIDITY: 
Tk=Te*#273.15 


Ki=2.28247@8E-1+Tk9(S.S989286E-4+Tke3 .S5714286E-7) 
K2=2 .SU46515E-l+1k (2 3978S / le -Seike (lees o7 lea 
K3=3.2138932E-1-Tke( 3.00428S57E-4-Tkel .408S57148-6 ) 


AZ=(K3-2¢K2+K1 1/200 

Al=(K2-K1-90Q*#A2 1/1 

AQ@=K1-40*eA1L-160G0*A2 

KmAQ+Egr*(AL+EgreAd ) 

RETURN K 

FNENO 

DEF FNTanh( xX) 

P@EXP(X ) 

Q=1/P 

Tanh=(P-Q)/(P+Q) 

RETURN Tanh 

FNENDO 

COEF FNTvsv(V) 

COM 7Ce7 ebay) 

T=aC(Q) 

FOR [=1 TO 7 
—«TaT4+C (Cl eV 

NEXT I 

RETURN T 

FNENO: 

OEF FN8eta(T) 

Rop=FNRho( T+. 1) 

Rom2FNRho( T-. 1) 


Betaz-2/(Roo+Rom)*(Rop-Rom)/.2 


RETURN 8eta 
FNENO 
COEF FNPsat(Tc) 


@ TO 80 deg F CURVE FIT OF Psat 


TFL. BeTc+32 


PaeS .9455S5+TFe(  15352082+T ft 1. 4ee0955E—s+1 f *o 6 le0s7 te one 


Py=Pa-14.7 

IF Pg-@ THEN 
Psat=Pg 

ELSE 
Psat=Pge29.92/14.7 

ENO IF 

RETURN Psat 

FNENO 

SU8 Purge 

BEEP 


INPUT “ENTER FILE NAME TO SE DELETED” Files 


PURGE Files 


-2@ TO 20 DEG C 


+=FS1G,-=1n Hg 


GOTO 2919 

SUBEND 

SUB Wilson 

COM /Will/ Doa(4),014a(4),Kma( 4) [act ,Oroot(4),La(Z) 
COM /Wil2/ Delta,Isat,Nsets,Hod,C1a(4),Alpaa(3) 

DIM Emf(17),TC17),Xa(100),Yal 109) ,Pe(4),Tp(4) 

BEEP 

INPUT “PLEASE RE-ENTER NAME OF FILE” ,Filed 

ASSIGN @File TO Files 


INPUT *@=SMOOTH COPPER TUBE, L#=CUNI FIN, c=KORODENSE, 3=FIN TITAN”, [tube 


BEEP 


INPUT “@=NOQ INSERT, != TWISTED TAPE INSERT, 2*HITRANI, 32HITRAN2Z*, Insert 


BEES 

INPUT “GIVE A NAME FOR XY FILE" ,Xy3 
CREATE BOAT Xy$,20 

ASSIGN @Xy TO XyS 

ey as a 

Do=Doal( [tube ) 

Dr.=Oi1a(I[ tube) 

Km@kKma( IT tube) 

Ax=PT*#0i*2/4 ! Cross-sectional area 
Ao=PI *Doel 

Rm=Qo*L0G(Do/Di )/(2*Km) 


Initial values 

Tf=Tsat 

Alpa=.6SS 

IF Insert=@ THEN 

Ciwl.O, 

ENO IF 

IF Insert#=l THEN 

Ci12#5.172 

ENO IF 

IF Insert=2 THEN 

Gree cb 

ENO IF 

IF Insert=3 THEN 

Ci=.063 

ENO IF 

6=9.31 

Ibeg=2 

lend=0 '!CHANGE TO 4, IF FIVE TUBES IN BUNDLE 

IF [Lact”1@ THEN 
[Ibeg2lact 
Tend=lact 

ENO IF 


FOR [=Ilbeg TO lend 
5420 
S720 
S$«5s=0 
Sxy 20 
Jset 20 
ASSIGN @File TO Files 
ENTER @File:Otges ltube,Egrat ,.Odl,Od2,0d3,0d4 ,OdsS 
ENTER @File:Pct(*),Emf(*), Tote) 
FOR J=@ TO 17 
TC J @FENTvsvi(Emf( J )) 
NEXT J 
Tvap2(7T(@)4T( 194T(2) 9043 


hss, 


3162 Tliag2(T(3)+T(4))/2 


S15! Tligt=T(S) 
3165S IF Isat=@ THEN 
mpl ya! Tsat#Tligl 
By ELSE 
3180 Tsat=Tvap 
3188 ENO IF 
3199 |! Grad=FNGrad( T(I+6)) 
oral Touts(T(2e1+10)+T( 2el+ll))/2 
Sys hs Delt=Tout-T(I+6) 
3200 Tavg=I( [+6 )+0elte.S 
3201 ! IF I#lend THEN 
3202 ! Oelt=T(10)-T(9) 
3204 | Tavga(T(10)-T(9))/2 
3205 |! ENO IF 
3206 ! 
S210! Water/Ethylene Glycol Mixture Properties 
3215 Rhoeg#FNRhoeg Tavg,Egrat ) 
3220 Nueg#FNNueg< Tavg,Egrat ) 
52a5 Mueg=Nueg*Rhoeg 
3230 Cpeg=FNCpeg! Tavg,Egrat ) 
Z2o5 Keg=FNKeg( Tavg,Egrat ) 
3248 Preg=Cpeg*Mueg/Keg 
3245 | 
S250 Mdot#FNFmecal(I,T(I+6),Pe(I)) 
5295 Veg=Mdot/(Rhoeg*Ax ) 
3260 Reeg=Veg*01/Nueg 
3265 Ress=4eMdot/(Mueg*(Pl1*01-490elta)) 
3278 Qdot=Mdot *Cpeg*Delt 
S215 Qdp=Qdot/Ao 
3280 Lmtd=Delt/LOG( (Tsat-T(I1+6))/( Tsat-T(I1+6 )-Delt)) 
S28s Vo=Qdp/Lmtd 
286 IF Insert2=@ THEN 
Soa Bbl=.0668*(01/L )*Reeg*Preg 
Saco Bbo2=1+.049((01/L )*Reeg*Preg)* .6666 
ves Omega=3.66+(8b1/Bb2) 
3290 ENO IF 
BAS h | IF Insert! THEN é 
Ba: ae: Omega2(1+S5.484E-3*Preg* .79(Res/Hod)*1.25)*.S5 
3294 ENO IF 
e225 IF Insert=2 THEN 
SESE Omega?(Reeg*.65)*(Preg” .46) 
5237, ENO IF 
gee IF Insert235 THEN 
yl S| Omega=(Reeg*.76)*(Preg*.46) 
3222 ENO IF 
3202 |! 
35@3 =| R-LL4/RIil3 Properties 
g505 Hf g=FNHfg(Tsat ) 
551 KF sFNK (CTF) 
gon Rhof=FNRho( TF) 
cee Muf=FNMu( TF ) 
Sos! 
3320 Fai KP" 3*Rhof*2*GeHfg/(Muf *Do*Qdp ) )*.33335 
LAO HomAlpaeF 
33549 Two=Tsat-Qdp/Ho 
$345 TRaTsat/st+leTwo/3 
3252 Y¥a(t/Uo-Rm) oF 
EoS5 X=Do*F/( Keg*Omega ) 
3362 | PRINT “OMEGA="; Omega: Fa". Fi "Xa" ixi "Ya"; y 
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Xa Jset )sx |  [NEFFICIENT (MODIFY LATER) 

Ya(Jset = 

Sx 25.4 +X 

Sy=Sy+Y 

Sxs2Sxs+X ox 

SxysSxytXeyY 

Jset=Jset+l 

IF Jset<Nsets THEN 3135 

ASSIGN @File TO » 

Slope=(NsetseSxy-Sx#Sy )/(NsetseSxs-Sx "2 ) 

Intcpt2(Sy-SlonpeeS« )/Nsets 

Ciczl/Slope 

Alpac=l/Intcpt 

Cerr=ABS((Ci1-Cic)/Cic) 

Aerr2ABS( (Alpac-Alpa)/Alpac ) 

IF Cerr>.@@l1 OR Aerr-.@@1 THEN 
Alpa=(AlpatAlpac)#.S 
CraCOrtCie oes 
PRINT °“CIC=";Cic: “ALPA=";Alpa 
GOTO 3120 

ENO IF 

BEEP 

BEEP 

PRINTER IS 1 

PRINT “CIC#=";Cic:; “ALPA=":Alpa 

Cia(Insert )2Cic 

Alpaa(I)#Alpac 

PRINTER IS 7@1 

FOR J=@ TO Nsets-l 

OUTPUT @Xyixat(J),Ya(J) 


NEXT J 
PRINTER [Sat 
NEXT I 
ASSIGN @Xy TO « 
SUBENO 


DEF FNFmcal(T,Pe) 
Mdot15=.0019*#Pc-.0098 
“Mdot@=.Q@Q02@6«Pc-.G0414 
Mdot24=.8@@225*Pc-.8@Q@329 

IF T<@ THEN 

Mdt=((T+#15)*(Mdot@-MdotlS)/1S )+MdotlS 

ELSE 

Mdt=( Te(MdotZ4—-Mdot@)/24 )+Mdot@ 

ENO IF 

RETURN Mdt 

FNENO 

DEF FNCoshn(X ) 

PsEXP(X’) 

QaEKP( =X ) 

Cosh=.S5e(P+Q) 

RETURN Cosh 

FNENO 


re 
mn 
oa) 
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